Universita degli Studi di Pisa 
Facolta di Scienze Matematiche, Fisiche e Naturali 
Corso di Laurea in Fisica 



Violazione di CP e rivelazione dei decadimenti 
del K°: studio e simulazione del Trigger neutro 
nell'esperimento NA48. 



Tesi di laurea di 
Federico Calzolari 



Relatore: Prof. Giuseppe M. Pierazzini 



Anno Accademico 1994-95 



Indice 



1 



Indice 



1 Introduzione 4 

2 Violazione di CP e misure sperimentali 7 

2.1 La violazione di CP nei sistemi di particelle K° 7 

2.2 La misura di R 13 

2.3 Misure di ^ 13 

2.4 Violazione di CP in altri sistemi fisici 14 

3 L'apparato sperimentale NA48 15 

3.1 I fasci di K s e K L 17 

3.2 Le anticoincidenze e i contatori di veto 18 

3.3 II calorimetro elettromagnetico 18 

3.4 L'odoscopio neutro 20 

3.5 Lo spettrometro magnetico 20 

3.6 II calorimetro adronico 20 

3.7 L'odoscopio carico 21 

3.8 II contatore muonico 21 

3.9 Trigger 21 

3.10 Una prima fase di trigger 22 

3.11 II Trigger carico 22 

3.12 II Trigger neutro 22 

3.13 II Trigger Supervisor 23 

4 Simulazione Montecarlo 25 

4.1 Decadimenti del K e rivelatori: Montecarlo 25 

4.2 Simulazione fisica 26 

4.2.1 I fasci K s e K L 26 

4.2.2 I modi di decadimento dei K 27 

4.3 I rivelatori 28 



Indice 2 



4.3.1 II sistema di etichettatura dei K§ 28 

4.3.2 Le anticoincidenze e i contatori di veto 28 

4.3.3 Lo spettrometro magnetico 28 

4.3.4 L'odoscopio carico 28 

4.3.5 II calorimetro adronico 29 

4.3.6 II calorimetro elettromagnetico 29 

4.3.7 L'odoscopio neutro 30 

4.3.8 II contatore di veto per i muoni 30 

4.4 Sciami elettromagnetici 30 

4.5 Trigger 31 

5 Decadimenti neutri di K$ e Kl 32 

5.1 Analisi relativistica del decadimento 32 

5.1.1 Massa invariante 32 

5.1.2 Distribuzione di gamma e pioni 33 

5.2 Analisi off-line 35 

6 Simulazione del Trigger neutro 36 

6.1 II rumore 36 

6.2 Amplificatori 36 

6.3 Shaper 37 

6.4 II clock 39 

6.5 Struttura pipelined 39 

6.6 I sommatori analogici 40 

6.7 I digitalizzatori 40 

6.8 I filtri 41 

6.9 I sommatori in righe e colonne 42 

6.10 II Peak-Sum 42 

6.11 II Peak-Finder 45 

6.12 La temporizzazione dell'evento 45 

6.13 Look-up table 46 

6.14 II Trigger Supervisor 47 

7 Test e risultati 48 

7.1 Campionamento del segnale ed errori di arrotondamento 48 

7.2 Distribuzione temporale degli eventi 49 

7.3 Analisi e acquisizione dei dati 50 

7.4 Energia 51 

7.5 Numero di gamma 53 



Indice 3 



7.6 Baricentro dell'energia 55 

7.7 Vertice di decadimento della particella 57 

7.8 Efficienza del Trigger neutro 59 

7.8.1 Variazione di soglia 64 

7.8.2 Variazione di raggruppamento delle celle 64 

7.9 Eventi accidentali 64 

Bibliografia 68 



Capitolo 1. Introduzione 



4 



Capitolo 1 



Introduzione 



Nell'ambito dell'esperimento NA48 in corso al CERN Super Proton Syncrothon (SPS) 
finalizzato a misurare la violazione diretta di CP in sistemi di particelle K° attraverso il 
doppio rapporto 



con una precisione di -Re(^) superiore a 2 x 10~ 4 [1], il lavoro di tesi si propone di effettuare 
una simulazione per la rivelazione dei decadimenti neutri dei K° ed una successiva analisi 
dei dati raffrontando l'evento flsico generato con metodo Montecarlo e quello ricostruito. 
Per la misura di - occorre isolare tra i vari modi di decadimento dei K quelli Kg — > 7r°7r°, 
K s — > 7r + 7r~, e quelli K L — > 7r°7r°, K L — > 7r + 7r~. 

Due fasci quasi collineari di K s e K L sono creati in modo tale da avere una comune 
regione di decadimento; contatori per l'identificazione dei K s e anticoincidenze per la 
soppressione dei decadimenti al di fuori della regione fiduciale sono disposti lungo la 
traiettoria dei fasci. Oltre la zona di decadimento e perpendicolarmente ai fasci si trovano 
il calorimetro elettromagnetico a Krypton liquido per la rivelazione dei fotoni, prodotti del 
decadimento di K° in 27r°(-^ 47) e 3n°(— > 67), e il calorimetro adronico per la rivelazione 
dei decadimenti in 7r + 7r~ . 

II problema principale per la rivelazione dei decadimenti K° — > 2ir riguarda la eli- 
minazione del fondo nei decadimenti Kl — > 2tt : questo e' costituito essenzialmente 
dai decadimenti — > TT^e^u, — > ti^^u, Kl — > 7r + 7r _ 7r° per il modo carico, e dal 
decadimento — > 7r 7r°7r° per il modo neutro. Per ridurre il fondo nei decadimenti carichi 
sono stati progettati uno spettrometro magnetico e un contatore di veto per i muoni. 
Compito del Trigger neutro e' invece discriminare i decadimenti K s — > 27r° e K L — >■ 2ir° 
da Kl — > 37r° (questi ultimi costituiscono piu' del 99% dei decadimenti neutri dei Kl), e 
arrivare alia ricostruzione di numero di gamma incidenti sul calorimetro elettromagnetico, 



[1.1] 
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energia, baricentro dell'energia e vertice di decadimento della particella K. Di qui la 
necessita' di rimuovere quanti piu' possibile decadimenti Kl — > 37r°, compresi quelli con 
soli 4 gamma sul calorimetro elettromagnetico, cosi' da ridurre il numero di dati da 
registrare per una successiva analisi off-line. Tutto l'apparato di Trigger neutro lavora 
separatamente sulle due proiezioni x e y ottenute dalla somma in righe e colonne delle 
128 x 128 celle del calorimetro elettromagnetico. 

Una vasta collaborazione CERN ha prodotto dati relativi a decadimenti K$ — > 27r°, 
Kl — > 2tt°, Kl — > 3tt° generati con metodo Montecarlo. II lavoro di tesi parte da una 
analisi di tali dati, cui e' stato necessario apportate alcune modifiche, e procede con: i) una 
simulazione completa dei dispositivi per il trattamento dei segnali provenienti dalle singole 
celle del calorimetro elettromagnetico e per la digitalizzazione dell'informazione, ii) una 
simulazione dell'analisi effettuata dal Trigger neutro, iii) una ricostruzione dei parametri 
dell'evento fisico in ingresso. La simulazione dei vari componenti e' stata effettuata in 
Fortran, al fine di ottimizzare i tempi necessari per l'analisi numerica di un evento e degli 
algoritmi impiegati per la ricostruzione dell'evento stesso. 

Attraverso una accurata analisi di decadimenti K s — > 2n°, K L — > 2ir°, K L — > 3n°, 
sono state ottenute le distribuzioni di energia, angolo di apertura e punto di impatto sul 
calorimetro elettromagnetico dei gamma prodotti dal decadimento dei 7r°, e, attraverso la 
ricostruzione della massa invariante dei 7r° a partire dall'energia depositata sulle singole 
celle, la traiettoria virtuale dei n°. 

Da una successiva simulazione della ricostruzione dell'evento fisico da parte della catena 
elettronica progettata in pipeline, e' stato possibile ottimizzare le soglie e gli intervalli 
di accettazione delle varie quantita' fisiche in gioco, in corrispondenza di diverse confi- 
gurazioni possibili dell'elettronica stessa. E' stata infatti testata l'efficienza del Trigger 
neutro e degli algoritmi impiegati in condizioni di differenti metodi di compattazione 
dell'informazione; tale compattazione si e' resa necessaria, nonostante una inevitabile 
perdita in precisione, per una analisi da parte dell'elettronica sui valori di rate massimi di 
funzionamento, previsti intorno a 40 MHz. Dalla conoscenza dell'evento fisico reale in in- 
gresso all'apparato, e dalla sua ricostruzione via hardware, sono state calcolate l'efficienza 
e la precisione del Trigger neutro, e raffrontate con le aspettative teoriche. 
Particolare attenzione e' stata riservata all'analisi del comportamento del Trigger neutro 
in condizioni reali di funzionamento, ovvero con possibili sovrapposizioni spaziali e tem- 
porali, parziali o totali, di eventi, e in presenza di eventi accidentali, siano essi dovuti a 
sorgenti esterne, o a modi di decadimento diversi da quelli analizzati. 

Per quanto riguarda la simulazione dell'apparato di Trigger neutro sono stati creati 
dei moduli Fortran che descrivono il comportamento dei vari blocchi: amplificatori e 
shaper collegati alle singole celle del calorimetro elettromagnetico per formare il segnale 
da analizzare, sommatori in blocchetti di celle per la compattazione dell'informazione, 
digitalizzatori (fast ADC), filtri spaziali e temporali per l'abilitazione alia somma in righe 
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e colonne dei singoli blocchetti, Peak-Sum e Peak-Finder progettati all'INFN sezione di 
Pisa per il calcolo nelle due proiezioni x e y del numero dei picchi, della temporizzazione 
fine dei picchi stessi, e dei momenti dell'energia di ordine 0,1,2 rispetto alia posizione 
proiettata, necessari alia Look-up- Table per la ricostruzione di energia totale depositata 
sul calorimetro, baricentro dell'energia, vertice di decadimento del K°. 

E' stata dedicata grande attenzione alia verosimiglianza con la realta' fisica, nonche' 
alia velocizzazione del sistema per poter simulare ed analizzare una quantita' di dati sum- 
cienti a garantire una valida statistica dei valori inerenti i tre diversi modi di decadimento 
K s — > 2n°, K L — > 2tt°, K l — > 37r°, in assenza ed in presenza di eventi accidentali. 
I risultati di una analisi completa, pur in accordo con le previsioni teoriche, suggeriscono 
piccole modifiche al settaggio di alcuni parametri per una migliore emcienza nella acqui- 
sizione dati. 
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Capitolo 2 

Violazione di CP e misure 
sperimentali 

2.1 La violazione di CP nei sistemi di particelle K 



Dalla scoperta della violazione di CP in un sistema di mesoni K neutri, nel 1964, da 
parte di J.H.Christenson, J.Cronin, V. Fitch e R.Turlay [2], svariati esperimenti si sono 
susseguiti al fine di una migliore comprensione della portata e deH'origine del fenomeno. 
Fino ad allora si era creduto che tutti i processi elementari fossero invarianti sotto la 
applicazione di ciascuno dei tre operatori C, P, T separatamente. Tuttavia, gia' nel 
1957, una serie di esperimenti mise in evidenza la violazione di C e P separatamente nei 
decadimenti deboli dei nuclei [3] e delle particelle 7r [4] e fj, [5]. 

II sistema dei mesoni K° e' finora stato l'unico in cui si sia potuta misurare con notevole 
precisone la violazione di CP: nel decadimento dei Kl a lunga vita media in una coppia 
di 7r neutri o carichi, nell'asimmetria di carica nei decadimenti semileptonici dei Kl . 

In un sistema K° -K° , ne' \K°) , ne' \K°) sono autostati di CP. Supponendo inizial- 
mente che le interazioni deboli siano assenti, gli stati \K°) e \K°) sono autostati della 
Hamiltoniana H delle interazioni forti piu' interazioni elettromagnetiche 



[2.1] 



Hq = Hf + H ( 



em 



[2.2] H \K°) = M \K°) 

dove M e' la massa del K°. 

Introducendo gli operatori coniugazione di carica C e parita' P, definiamo: 



[2.3] 



\K°) = CP\K°) 
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Dal momento die H commuta con CP 
[2.4] H \W) = M \W) 



particella 


massa [CreV] 


K° 


0.497 
0.497 



Tab. 2. 1.1: Massa invariante dei K 

Inserendo neH'Hamiltoniana H le interazioni deboli (Hd), responsabili dei decadimenti 
di \K°) e \K°) , 

[2.5] H — Hf + H em + E d 

Noto che H d non conserva ne' C, ne' P separatamente, e' comunque lecito supporre 
inizialmente che Hd conservi il prodotto CP. 

Gli autostati simultanei di H e CP possono essere costruiti attraverso combinazioni di 
\K°) e \K°) : 

[2.6] 



1 



[2.7] 



\Ko 



lo stato \Ki) con CP = +1, lo stato \K 2 ) con CP = — 1. \Ki) e |i^ 2 ) sono rispettivamente 
gli stati a vita media breve e lunga. 

Considerando i decadimenti del K° in 2n o 37r , poiche' lo stato 27r ha CP = +1, il 
decadimento in In e' consentito solo per la particella K\ , e proibito invece per la particella 
K 2 , che quindi puo' decadere esclusivamente in 37r . Poiche' inoltre il volume nello spazio 
delle fasi associato al decadimento di un K° in 2n e' molto maggiore di quello associato 
ad un decadimento in 37r , la vita media del K\ e' inferiore a quella del K 2 : 



[2.1 



particella 


vita media 


s 




^2 


0.89- 10~ 10 
5.18 • 10" 08 



Tab.2.1.2: Vita media dei K 



Nel 1964 J.H.Christenson, J.Cronin, V. Fitch e R.Turlay [2], utilizzando un apposito 
spettrometro a due bracci per la rivelazione dei decadimenti a 2n munito di rigeneratore, 
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scoprirono che anche il K° a lunga vita media (K 2 ) puo' decadere in 2ti , e piu' precisa- 
mente in 7r + 7r~ con un branching ratio di 0.002, in 7r°7r° con un branching ratio di 0.001, 
violando percio' la conservazione di CP. 

Nel 1965 un esperimento di V. Fitch et al. [6] dimostro', attraverso l'interferenza co- 
erente tra le ampiezze di decadimento di K L e K s in 7r + 7r~ nel vuoto e in presenza di 
un rigeneratore, che i 2n carichi prodotti dal decadimento dei K L sono identici a quelli 
prodotti dal decadimento di Ks , eliminando quindi tutti i tentativi di spiegazione alter- 
nativa che avevano seguito il primo esperimento del 1964, quale tra gli altri la presenza 
di qualche altra particella a energia molto bassa. 

Gli stati con vita media e massa definita (\Ki) e \K 2 ) nel caso di conservazione di CP) 
saranno percio' i nuovi stati definiti \Ks) e \Kl) , rispettivamente a corta vita media e a 
lunga vita media. Utilizzando la notazione introdotta da Wu e Yang [7], gli stati \K$) e 
\Kl) risultano dalla combinazione di stati \Ki) e \K 2 ) : 



[2.9] 



\K s ) = ^=(\K 1 ) + e\K 2 )) 



[2.10] 



\K L ) = -j={\K 2 )-e\K 1 )) 



dove e e' un parametro complesso (|e| m 2 • 10 -3 ) che tiene conto della violazione di CP. 

La differenza di massa Mk l — Mk s = 3.521 • 10~ 15 GeV e' stata ottenuta attraverso 
esperimenti di decadimenti in 2n di miscele di K s e K L . 

Data una generica sovrapposizione degli stati \K°) e \K°) , l'evoluzione temporale del 
nuovo stato risulta data da: 



[2.11] 

[2.12] 
[2.13] 



$(t) = a(t)\K°) + b(t)\K°) = 



d_ 
di 



a(t) 
b(t) 



= l -x 

i 



a(t) 
b(t) 



L'operatore X,Hermitiano, puo' essere scomposto in: 

r 



[2.14] 



X = M-i- 



a(t) 
b(t) 



dove M e T sono rispettivamente la matrice di massa e la matrice di decadimento. 
Gli elementi della matrice X sono: 



[2.15] 



{K°\H\K°) 
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[2.16] (K \H\K°) 
[2.17] (K°\H\K°) 
[2.18] (K°\H\K°) 
L'invarianza sotto CPT impone die: 

[2.19] X n =X 22 

quindi 

[2.20] M u = m 22 r n = r 22 

mentre per il fatto che l'operatore e' Hermitiano: 



[2.21] 



X\2 — X 2 i 



quindi 

[2.22] M 12 = M* 21 T 12 = T* 21 

Le due matrici M e T risultano quindi essere: 



[2.23] M = 

[2.24] r = 



con M e T reali. 

L'equazione agli autovalori per l'operatore X ha come soluzioni: 



M 


M 12 


M* 2 


M 


r 


r i2 " 


1 12 


r 



[2.25] 
[2.26] 

relative agli autostati 
[2.27] 

[2.28] 



Ms — -T s = Xu + \/X 12 X 21 



Mr - ^Tr = X 



11 



X\ 2 X 2 i 



\Ks) 



\K T ) = 



1 



1 + |e| 
1 

1 + lei 



-AKJ + e\K 2 )) 



:(\K 2 ) + elK,)) 



dove e' stato introdotto il parametro e 



[2.29] 



V X\ 2 — \[X- 



21 



\JXy2 + \JX 2X 
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Sostituendo a i valori ricavati dall'equazione agli autovalori, si ottiene: 

^F 12 + ^M 12 

[ ■ 1 i(r s — r L ) — (m s — m l ) 

Utilizzando i dati sperimentali per Ms, Ml, Ts, IV, il valore del modulo di e risulta essere: 

[2.31] |e| = (2.3 ±0.02) • 10~ 3 

trascurando QT 12 , la fase di e e' determinata dal denominatore, ed e' pari a: 

[2.32] arg(e) = tg' ^ 1 ~ 5 S » 45° ± 1.5° 

II parametro e tiene conto della violazione indiretta di CP, attraverso il decadimento 
di una miscela di stati di opposta CP (\Ki) e \K 2 ) ). 

Nei decadimenti di K° in 2n lo stato finale puo' essere caratterizzato dal numero quan- 
tico di isospin 1 = 0, oppure 1 = 2, entrambe combinazioni simmetriche di due stati a 
isospin 1 = 1. Indicando con |0) e |2) gli stati di isospin corrispondenti rispettivamente a 
I — e / = 2 per i decadimenti in 2n , si ricavano quattro ampiezze di transizione: 

[2.33] (0\T\K s ) 

[2.34] (2\T\K S ) 
[2.35] (0\T\K L ) 
[2.36] @\T\K L ) 

dove T e' l'Hamiltoniana responsabile dei decadimenti che violano CP, riducibili a tre 
quantita' normalizzando rispetto all'ampiezza (0\T\K s ): 

(0\T\K L ) 



[2.37] e = 

[2.38] el = 

[2.39] to = 



(0\T\K s ) 

(2\T\Kl) 
(Q\T\K S ) 

Q\T\Ks) 



(0\T\K s ) 

Le quantita' osservabili, disponibili da esperimenti, sono: 

7i + 7i-\T\K L ) e + e/ 



[2.40] ^ = 



(n+n-\T\K s ) 1 + 



V2 
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[2.4!] „„ ^™ f " 2 " 



(7i°7i \T\K s ) l-y/2u> 

r 2 42l r (K L -> TC+e-u) - (Kl ~> 7T£+V) 

[ ' J L (K L ^ 7T+£-u) + (K L -> 7T-£+Z/) 

dove £ e' un leptone: elettrone(positrone) o muone z^. 

Definendo gli stati a due pioni 7r in funzione degli autostati dell'operatore isospin: 

[2.43] 

[2.44] 
poiche' w < 1 : 

[2.45] 7]+ _ = e + el 

[2.46] 7700 = e - 2e/ 

[2.47] 5 L = 23fc 

La grandezza fisica misurabile 5l e' riferita ad una asimmetria di carica nei decadimenti 
semileptonici if L — > tt^Pu, ed evidenzia anch'essa la violazione di CP nei sistemi di K° . 
Introducendo le ampiezze di transizione di K° e K° nei canali a isospin I = e / = 2: 




[2.48] 


(0|T|fsT ) = A e i<5 ° 


[2.49] 


(2\T\K°) = A* Q e iSo 


[2.50] 


(0\T\K°) = A 2 e iS2 


[2.51] 


(2\T\K°) = A* 2 e i52 



il parametro el risulta essere 

[2.52] el = ^Ae^- 5 ^ 

ed e' direttamente collegabile alia violazione diretta di CP, attraverso il decadimento 
diretto di uno stato a CP definita in uno stato finale a CP opposta. 

La scomposizione delle ampiezze di decadimento dei K° osservabili in e e el corrisponde 
ad una separazione degli effetti di violazione di CP a livello della matrice di massa M (e) 
e a livello della matrice di decadimento T (el). 

II doppio rapporto 

2 



[2.53] R 



V+- 



N(K L - nW) f N(K L ^n+n-) ^ ^ _ ^ 



N(K S -> TT°TT )' N(K S -> 7T + 7T-) V 

costituisce la misura cui e' finalizzato l'esperimento NA48 al CERN Super Proton Syn- 
chroton (SPS). 
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2.2 La misura di R 

II modulo e la fase di r/ + _ sono stati misurati in maniera precisa attraverso fasci di K s e 
Kl preparati attraverso un rigeneratore: 

[2.54] r] + _ = [(2.274 ± 0.0022) • 10~ 3 ] exp[i(0.778 ± 0.021)] 

Per quanto riguarda la misura di modulo e fase di i] Q0 , questa risulta piu' complessa da 
determinare, e comporta errori sistematici maggiori che nella misura di questo perche' 
dopo il decadimento del K in 7r°7r° , il ir° decade con vita media brevissima (dell'ordine 
di 10~ 16 s) in due fotoni, per cui ne risulta un evento a quattro fotoni, energia e posizione 
di ciascuno dei quali devono essere misurate con estrema precisione. 



particella 


vita media 


S 




massa [CreV] 


7T + 
7T~ 
7T° 


1Q -08 
1Q -08 
1Q -16 


0.139 
0.139 
0.135 



Tab. 2. 2.1: Vita media e massa dei n 



2.3 Misure di eJ - 

e 

La prima prova della violazione diretta di CP e' stata annunciata della collaborazione 
NA31 del CERN di Ginevra nel 1988 [8], con un valore 

[2.55] SR(^) = (2.30 ±0.65) • 10~ 3 

Nel 1990 l'esperimento E731 al Fermilab di Chicago [9] ha riportato il valore 

[2.56] SR(^) = (0.74 ±0.59) • 10~ 3 

L'attuale migliore stima di ^ dovrebbe quindi risultare dalla media, pesata rispetto al 
numero di decadimenti analizzati, di queste due precedenti misurazioni: 

[2.57] £(-) = (2.10 ±0.90) • 10~ 3 

Nell'esperimento NA31, concettualmente simile a NA48, salvo una minore precisione, 
la precisione statistica e' stata limitata dal numero relativamente basso di decadimenti 
K L -> 7T 7r° . 

Altri esperimenti sono stati effettuati negli ultimi anni, senza pero' fornire misure con 
precisione tale da poter essere confrontate con queste. 
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2.4 Violazione di CP in altri sistemi fisici 

Sistemi analoghi al sistema K° -K° sono D°-D° e B°-B°, cui e' possibile applicare lo 
stesso formalismo impiegato per i mesoni K neutri. Tali sistemi sono pero' costituiti da 
particelle molto piu' pesanti, con vite medie molto piu' brevi, e con un vasto numero di 
possibili stati finali di decadimento. 

Nonostante in questi sistemi la violazione di CP a livello della matrice di massa M sia 
maggiore che nel caso del sistema K° -K° , risulta molto piu' difficile da misurare. Nel 
sistema K° -K° infatti, essendo la massa del K° poco piu' grande della massa totale dei 
3n , il che implica una lunga vita media del Kl , esiste una grande differenza tra la vita 
media delle particelle Ks e Kl , consentendo quindi una buona identificazione, all'atto 
dell'esperimento, delle particelle stesse. 
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Capitolo 3 

L'apparato sperimentale NA48 



L'esperimento NA48 [1], attualmente in corso al CERN Super Proton Synchroton (SPS), 
consta di un apparato di rivelazione per la misura della violazione diretta di CP nei 
decadimenti dei K neutri, attraverso il doppio rapporto 

con una precisione di Re(^) superiore a 2 x 10~ 4 . 

La misura di - e' ottenuta attraverso un raffronto dei rates relativi di decadimento di 

e 

Ks e K L in n + 7i~ e 7r°7r° . 



Fig.3.0.1: Apparato sperimentale NA48 
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Per una misura precisa occorre che le efficienze di rivelazione per il modo carico e per il 
modo neutro risultino il piu' possibile uguali; questo e' ottenuto attraverso l'utilizzo degli 
stessi contatori per la rivelazione dei decadimenti in 7i + ir~ e n 7i° , siano essi derivati da 
K s , sia da K L . 

In tabella sono riportati i principali modi di decadimento dei Ks e Kl [12]: 



particella 


decadimento 


branching ratio 


K s 


+_- 

7T 7T 


n pop 




7T 7T 


0.314 




7r + 7T _ 7 


0.002 




7T 7T 71 


0.005 • 10~ 2 






0.004 • 10~ 2 




e + e~ 


0.001 • KT 2 




77 


0.002 • 10~ 3 






0.003 • 10" 4 


K L 




0.386 






0.270 




ttVtt 


0.217 






0.124 






0.013 






0.002 




ttV 


0.009 • KT 1 




77 


0.005- 10- 1 




7T°77 


0.002- 10- 1 




7r°7r ± e T v 


0.006 • 10- 2 




7T + 7T _ 7 


0.004 • 10~ 2 



Tab.3.0.1: Modi di decadimento dei K neutri 



II problema principale riguardante la rivelazione dei decadimenti dei K° consiste nella 
eliminazione del fondo per i decadimenti K L — > 2n . Questo e' costituito essenzialmente 
dai decadimenti K L — > 7r ± e =F z/, K L — > 7r ± / u =F z/, K l — > 7r + 7r _ 7r° per il modo carico, e dal 
decadimento K^ — > 7r 7r 7r per il modo neutro. 

Uno spettrometro magnetico e' progettato per ridurre il fondo nei decadimenti carichi, 
ed un calorimetro elettromagnetico a Krypton liquido per ridurre il fondo nei decadimenti 
neutri. 
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3.1 I fasci di K s e K L 

Disponendo di due fasci quasi collineari di K s e K L aventi una comune regione di decadi- 
mento, all'interno di una zona di 120 m di lunghezza e 2.5 m di diametro, mantenuta 
sotto vuoto spinto (Vacuum tank), e' possibile confrontare direttamente i rate relativi dei 
decadimenti di Ks e Kl in 7i + n~ e 7r°7r° . L'angolo di convergenza dei due fasci di Ks e 
Kl e' inferiore alia divergenza del singolo fascio, al fine di eliminare tutte le possibili fonti 
di differenza sistematiche nella rivelazione dei modi di decadimento di Ks e Kl . 

Un fascio di protoni primari con intensita' di 10 12 protoni / s ed energia di 450 GeV 
viene utilizzato per produrre il fascio di K L dall'impatto su una targhetta K L . I K L cosi' 
creati passano attraverso un primo collimatore ed un primo contatore, mentre i protoni e 
le altre particelle cariche, frutto dell'interazione dei protoni con la targhetta Kl , vengono 
deviati dalla direzione del fascio K L . 

Dopo un secondo collimatore, una piccola frazione dei protoni primari di intensita' pari a 3- 
10 07 protoni/ s viene nuovamente devialta verso il fascio dei K L e, attraverso un deflettore 
magnetico, viene convogliata lungo la beam-line costituita dal fascio Kl . Un ultimo 
deflettore devia nuovamente il fascio di protoni e lo indirizza, attraverso un contatore, 
verso la targhetta Ks per la generazione dei Ks ■ 

La differenza angolare nella traiettoria del fascio di Ks rispetto al fascio di Kl e' di 
0.6 mrad. 

In questo modo il numero di decadimenti Ks — > 2ir risulta essere confrontabile col numero 
di decadimenti Kl — * 2tt all'interno della comune regione di decadimento. 

La distinzione tra decadimenti di Ks e decadimenti di Kl e' affidata ad un sistema di 
tagging-counter per l'identificazione e l'etichettatura dei protoni impiegati per la genera- 
zione dei K s ■ Questa identificazione viene effettuata misurando la differenza temporale 
tra il passaggio dei protoni nei contatori posti a monte della targhetta K s e il tempo 
dell'evento nel rivelatore. Gli eventi per cui questa differenza temporale risulta essere 
all'interno di un intervallo di accettazione vengono etichettati come Ks , mentre gli altri 
vengono identificati come Kl ■ 

Qualunque inefficienza nel sistema di conteggio dei K s implica quindi che un decadi- 
mento Ks venga identificato come un decadimento Kl ; d'altra parte, ogni conteggio 
accidentale nel contatore Ks fa si' che un decadimento Kl venga identificato come un 
decadimento Ks ■ 

Importante e' mantenere basse queste probabilita' di transizione, e questo avviene at- 
traverso la scelta di una finestra temporale di accettazione di 5 ns. Cosi' facendo l'errore 
nell'identificazione di Ks e Kl e' ridotto a circa il 5%. 

II rate totale degli eventi nel rivelatore e' previsto intorno a 10 6 Hz. 
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3.2 Le anticoincidenze e i contatori di veto 

I decadimenti K L — > 7r°7r 7r° avvengono con un rate 200 volte maggiore dei decadimenti 
K L — > 7r°7T utili alia misura di R. E' percio' essenziale che una fase di pre-trigger fornisca 
un veto per la maggior parte di tale fondo, al fine di bloccare l'acquisizione dati per questi 
eventi. Parte di questo compito e' affidato ad una serie di anticoincidenze, costituite 
da un doppio strato di scinitillatore preceduto e intercalato da strati di Piombo per la 
rivelazione dei gamma, e posizionate in maniera tale da coprire un ampio angolo solido 
esterno al rivelatore. II cono di accettanza dell'apparato risulta cosi' delimitato da una 
serie di sette anticoincidenze disposte ad anello intorno alia beam-line, in grado di rivelare 
i gamma con traiettorie esterne ai rivelatori. 

L'efncienza di tali anticoincidenze per i decadimenti Kl — > 7r°7r 7r° e' pari a 0.75 per 
un'energia del K L di 100 GeV, e sale fino a 1.00 per K L ad energia inferiore a 50 GeV, 
che risultano quindi avere un angolo di decadimento dei n° , quindi dei gamma, maggiore. 
In questo modo viene soppressa parte parte del fondo per idecadimenti neutri in 7r°7r 7r° . 

Per il modo carico e' necessario sopprimere il conteggio per i decadimenti Kl — > 7i ± e T u, 
Kl n ± fi T i', Kl 7r + 7r~7r°. I decadimenti Kl — > 7r + 7r~7r° vengono anticoincisi at- 
traverso la rivelazione dei 7r + 7r~ sul calorimetro adronico e del 7r° sul calorimetro elet- 
tromagnetico; i decadimenti Kl — > TT ± e T u vengono anticoincisi dalla ricostruzione sul 
calorimetro adronico della massa invariante di un solo tt° ; i decadimenti K L — ► 7r ± /i T z/ 
vengono anticoincisi attraverso l'utilizzo di quattro larghi piani di scintillator! utilizzati 
come contatori di veto per i muoni fj, ± . 

3.3 II calorimetro elettromagnetico 

Per la rivelazione dei decadimenti neutri del K° e' necessario disporre di un rivelatore 
rispondente alle sequenti caratteristiche: 

• possibilita' di registrare eventi a piu' fotoni ad un rate di 1 MHz 

• risoluzione in energia a(E) = 0.035^ 

• risoluzione spaziale inferiore a 1 mm 

• risoluzione temporale inferiore a 1 ns 

Per soddisfare tali richieste e' stato progettato un calorimetro elettromagnetico ottagonale 
di 256 cm di raggio ad alte efficienza e risoluzione. 
Aspetti peculiari di tale apparato di rivelazione sono: 

• segmentazione del rivelatore in una matrice di 128 x 128 celle di 2 cm x 2 cm ciascuna 
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• bassa capacita' dei singoli canali del rivelatore, per permettere una veloce lettura 
dell'energia depositata 

• utilizzo di Krypton liquido ad alta densita' in una struttura quasi omogenea 

Ogni singola cella ha una superficie utile di 2 x 2 cm 2 , ed e' costituita da un elettrodo 
centrale ad alta tensione inserito tra due elettrodi di massa distanti 1 cm da esso e condivisi 
da due celle adiacenti sul piano orizzontale. 

Gli elettrodi, di una lega di Rame e Berillio, sono divisi in strisce la cui altezza (2 cm) 
definisce la struttura verticale della singola cella. Tali strisce, lunghe 125 cm, sono state 
disposte a zig-zag con angoli di ±50 mrad dall'asse della beam-line, al fine di evitare che un 
gamma incidente sullo spessore (40 /im)della singola striscia di rame possa incanalarvisi 
senza interagire col mezzo rivelatore, e quindi non rilasciare energia in nessuna cella. 

Come mezzo ionizzabile viene impiegato Krypton liquido (LKr) per una profondita' 
pari a quella degli elettrodi di rivelazione (125 cm), corrispondente a circa 27 lunghezze 
di radiazione. 

All'interno della singola cella ogni drift-gap e' di 1 cm, corrispondente ad un drift-time 
di 2.8 fis per una tensione applicata di 5000 V [13]. 

Uno sciame elettromagnetico, prodotto dall'interazione di un gamma con la materia 
(Krypton liquido), che si sviluppa all'interno del calorimetro, vi deposita una certa quan- 
tita' di carica che induce un impulso di corrente in corrispondenza di ciascuna cella del 
rivelatore. 

Tale impulso, di forma triangolare, viene raccolto da amplificatori di carica, direttamente 
connessi ad ogni singola cella, e contenuti all'interno del criostato per il mantenimento 
del Krypton alio stato liquido. 

II segnale lascia poi il criostato ed e' filtrato attraverso uno shaping-amplifier, che provvede 
a formare un impulso in uscita di ampiezza proporzionale alia derivata del segnale in in- 
gresso, e caratterizzato da una durata di 110 ns, con larghezza a meta' altezza di 60 ns, 
ed undershoot del 3% della durata di 2.5 fis. La raccolta di un'energia pari a 1 GeV sulla 
singola cella fornisce una corrente con un'ampiezza massima di 3 fiA, ed una successiva 
ampiezza del segnale di shaper di 10 mV. L'energia massima depositabile sulla singola 
cella e' di circa 50 GeV. 

La risoluzione in energia e' data da [13]: 



[3.2] 




0.035 



) 2 + 0.005 2 



e il profilo dell'impulso di corrente in funzione del punto di impatto del gamma sulla 
singola cella, a parita' di energia depositata, risulta essere quasi piatto. 
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La risoluzione spaziale e' inferiore a 2 mm, ed e' data da: 



[3.3] 



4.2 

a x = i = + 0.6 mm 



y/E [GeV] 



[3.4] 




3.4 L' 



odoscopio neutro 



Oltre al calorimetro elettromagnetico, un odoscopio neutro a scintillazione ha il compito 
di rivelare il tempo esatto di arrivo dell'evento, cosi' da costituire una fase di pre-trigger 
per i decadimenti carichi. 

Tale odoscopio, diviso in due semipiani, richiede una coincidenza destra/sinistra per 
l'accettazione degli eventi. 

3.5 Lo spettrometro magnetico 

I momenti delle particelle cariche vengono misurati in uno spettrometro magnetico costi- 
tuito da due serie di quattro camere a drift, tra le quali e' interposto un dipolo magnetico 
centrale. 

Con un campo magnetico integrato equivalente ad una variazione di momento trasverso 
di 200 MeV I c, la risoluzione nella misura del momento di una particella e' limitata dallo 
scattering multiplo sul materiale costituente le camere a: 



Questo e' in accordo con la scelta di camere a drift aventi una risoluzione spaziale di 
100 fim. 

Ciascuna delle quattro camere a drift consiste di quattro coppie di piani di fill formanti 
due sistemi di assi cartesiani ruotati fra loro di 45° per una migliore precisione nella 
misura. Ciascun piano e' costituito da 240 fill distanti fra loro 1 cm. 

Le camere a drift lavorano con una miscela 70/30 di Argon-isoButano, e la velocita' di 
drift risulta di 5 cm/ {is. 

3.6 II calorimetro adronico 

Oltre lo spettrometro magnetico e' disposto un calorimetro adronico ad Argon liquido 
per la rivelazione della posizione dei n + e n~ che, avendo vita media piu' lunga dei tt° , 



[3.5] 



Ap 

— = 0.006 
V 
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arrivano pressoche' tutti a colpire i rivelatori senza decadere. 



particella 


vita media 


S 




7T + 
7T~ 
7T° 


10 -u« 

1Q -08 
1Q -16 



Tab.3.6.1: Vita media dei it 



Dal punto di impatto del it sul calorimetro adronico e dalle informazioni ricevute dalle 
camere a drift, e' possibile ricostruire l'esatta traiettoria del 7r ± , e, dal raggio di curvatura 
impresso dal campo magnetico B dello spettrometro magnetico, l'impulso della particella. 

3.7 L'odoscopio carico 

Una fase di pre-trigger per i decadimenti carichi viene effettuata da un odoscopio posto 
di fronte al calorimetro elettromagnetico. Tale odoscopio, diviso in quattro quadranti, 
richiede una coincidenza in due quadranti opposti per l'accettazione dell'evento 

3.8 II contatore muonico 

Uno scintillatore costituito da 3 piani paralleli, posto oltre il calorimetro elettromagnetico 
e il calorimetro adronico, ha il compito di rivelare i decadimenti di Ks e Kl in cui sia 
presente un mesone /x, nonche' la posizione del punto di impatto del fi stesso, al fine di 
porre un veto sull'accettazione di tali eventi. 

3.9 Trigger 

Mentre i Ks decadono quasi esclusivamente in -k + tt~ e 7r°7r° , per quanto riguarda i Kl , 
questi decadono in 2n soltanto nello 0.3 % dei casi, e risultano pertanto coperti dal fondo 
degli altri predominanti modi di decadimento. 

Un complesso sistema di trigger a piu' livelli [14] e' necessario ai fini dell'esperimento 
per la rivelazione dei decadimenti in 2n , siano essi dovuti al decadimento di un Ks , o 
di un K L . Tenuto conto dell'alta intensita' del fascio, al trigger e' richiesto di avere una 
buona efncienza a rates elevati, e di riuscire a eliminare una elevata percentuale del fondo. 

I candidati per i decadimenti in 7r 7r sono selezionati tra quegli eventi con una energia 
ricostruita tra 60 GeV e 180 GeV, con un vertice ricostruito tra 2.1 m e 48.9 m dalla 
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posizione dell'ultimo collimatore, e un numero di gamma sul calorimetro elettromagnetico 
compreso tra 3 e 4. 

I candidati per i decadimenti in 7r + 7r~ sono selezionati tra i vari eventi utilizzando gli 
stessi tagli in energia, vertice e numero di gamma adottati per i decadimenti neutri [8]. 

3.10 Una prima fase di trigger 

Una fase di pre-trigger viene effettuata, al fine di ridurre il rate degli evnti da analizzare 
di circa 10 volte, dal Trigger di livello 1. Elaborando in tempo rapido le informazioni 
provenienti da anticoincidenze, contatori di veto, odoscopio carico ed odoscopio neutro, 
tale trigger abilita o meno la successiva analisi da parte dei trigger di livello 2: il Trigger 
carico e il Trigger neutro. 

3.11 II Trigger carico 

II trigger carico, utilizzato per distinguere i decadimenti 7r + 7r~ dal fondo, utilizza, per la ri- 
costruzione dell'evento, le informazioni provenienti dagli odoscopi carichi, dal calorimetro 
adronico, e dallo spettrometro magnetico. 

Per i decadiment carichi e' richiesta nei rivelatori una coincidenza simmetrica rispetto alia 
beam-line. Dalle informazioni del calorimetro adronico e dello spettrometro magnetico, 
il trigger carico ricostruisce energia, massa invariante, momento trasverso e vertice di 
decadimento del K . 

3.12 II Trigger neutro 

Compito principale del Trigger neutro e' la rimozione del fondo dominante costituito dai 
decadimenti Kl — > 7r 7r°7r° nei decadimenti neutri. 

I 7T° decadono con vita media brevissima generando ciascuno due gamma, che vengono 
poi rivelati sul calorimetro elettromagnetico. I decadimenti Kl — > 7r 7r°7r° si presentano 
come eventi a 6 gamma; tali eventi possono essere rivelati sul calorimetro come decadi- 
menti a 6 gamma nei caso in cui tutti i fotoni raggiungano il rivelatore, come decadimenti 
con numero di gamma inferiore a 6 nei caso in cui alcuni fotoni seguano traiettorie esterne 
all'apparato di rivelazione o vengano persi in quanto poco energetici. 

I tagli principali che consentono una notevole riduzione del fondo Kl — > 7r 7r°7r° vengono 
effettuati sul numero di picchi, sull'energia totale rilasciata nei calorimetro elettromagne- 
tico, sul baricentro dell'energia, sul vertice di decadimento del K . Al Trigger neutro 
e' richiesto di lavorare in modo conservative rispetto al numero di eventi, ovvero senza 
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scartare eventuali eventi sporcati da un evento accidentale poco energetico, die potrebbero 
magari essere ricostruiti ugualmente bene attraverso una successiva analisi off-line degli 
eventi registrati. 

II Trigger neutro lavora non sulle informazioni provenienti da tutte le 13500 celle del 
calorimetro elettromagnetico, bensi' sulle proiezioni x e y ottenute dalla somma in righe 
e colonne, ulteriormente accorpate in gruppi di due, dell'energia depositata sulle singole 
celle. 

All'interno di ogni singola riga (colonna) le celle vengono sommate in blocchetti 8x2 (2x8) 
e quindi digitalizzate. I blocchetti cosi' creati vengono analizzati dal filtro temporale, il 
quale abilita alia somma in righe e colonne il singolo blocchetto, quello temporalmente 
seguente e i due precedenti, nel caso in cui il blocchetto in esame risulti maggiore di una 
soglia selezionabile dall'esterno. 

I vantaggi di una siffatta compattazione dell'informazione, nonostante una inevitabile 
perdita in precisione ed efficienza, consistono nella possibilita' da parte dell'apparato 
elettronico di un'analisi sui valori di rate massimi di funzionamento, previsti intorno a 
40 MHz. 

Per una maggiore velocita' di analisi, l'elettronica e' stata progettata in pipeline, con tempi 
di elaborazione e di transito dell'informazione nel singolo stadio ampiamente inferiori ai 
25 ns, al fine di garantire una stabilita' dei registri al momento della commutazione del 
clock. 

Dispositivi principali dell'apparato di Trigger neutro sono il Peak-Sum e il Peak-Finder, 
progettati all'INFN sezione di Pisa per il calcolo nelle due proiezioni x e y del numero di 
picchi in energia presenti sul calorimetro elettromagnetico, della temporizzazione fine dei 
picchi stessi, e dei momenti dell'energia di ordine 0,1,2 rispetto alia posizione proiettata, 
necessari ai moduli Look-up Table per la ricostruzione di energia, baricentro dell'energia, 
vertice di decadimento del K . 

3.13 II Trigger Supervisor 

Per gli evnti giudicati buoni dal trigger carico e dal Trigger neutro, una successiva e piu' 
accurata analisi, resa possibile dai tagli effettuati e quindi da un rate di eventi decisamente 
inferiore, mira ad una migliore e piu' precisa ricostruzione delle quantita' fisiche osservabili, 
al fine di attivare l'acquisizione dati su disco per gli eventi candidati ad essere frutto di 
decadimenti K$ — > In e Kl — > 2n . 

Oltre alle quantita' precedentemente calcolate, ovvero numero di picchi, energia, bari- 
centro dell'energia, vertice di decadimento del K , il trigger supervisor provvede ad una 
accurata ricostruzione della massa invariante della particella originaria, e puo' quindi porre 
nuovi tagli su questa quantita', riducendo cosi' ulteriormente il fondo dei decadimenti del 
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K a tre corpi. 

II Trigger Supervisor provvede anche ad un raffronto incrociato delle informazioni prove- 
nienti dal trigger carico e dal Trigger neutro per una analisi dei decadimenti misti (in 
particelle cariche e neutre insieme). 

Alia fine dell'analisi effettuata on-line dall'apparato sperimentale, il fondo previsto per 
i quattro modi di decadimento utili alia misura e' [8]: 



decadimento 



fondo residuo (%) 




0.03 
0.07 
0.63 
2.67 



Tab. 3. 13.1: Fondo residuo 



II fondo residuo maggiore risulta essere, come previsto, per il modo di decadimento 
K L -> ttV . 
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Capitolo 4 

Simulazione Montecarlo 



4.1 Decadimenti del K e rivelatori: Montecarlo 

Al fine di studiare i vari aspetti del nuovo esperimento NA48 EPSI inerente la violazione 
di CP nei decadimenti dei K neutri, e' stato creato da una vasta collaborazione CERN 
un programma di simulazione per generare in tempi relativamente brevi una vasta gamma 
di decadimenti con metodo Montecarlo. 

II programma, denominato appunto New Montecarlo (NMC) [11], genera eventi di decadi- 
mento di Kg e Kl partendo da una banca dati (particle bank), relativa alia particella ori- 
ginaria, dove sono contenute informazioni quali tipo di particella, quadriimpulso, vertice 
di decadimento, energia, vita media. 

NMC simula i processi di decadimento di tali particelle e la risposta dei rivelatori uti- 
lizzati per le misure. Attraverso banche dati contenenti informazioni relative a traiettorie 
e modi di decadiemnto delle varie particelle, vengono ricostruiti i percorsi delle particelle 
e dei loro prodotti di decadimento fino ai rivelatori, disposti lungo la traiettoria del fascio 
e perpendicolarmente ad esso. 

Dopo il decadimento dei 7r° in gamma nella camera a vuoto per il modo di decadi- 
mento neutro, viene simulata la rivelazione dei gamma da parte del calorimetro elettro- 
magnetico a Krypton liquido attraverso l'utilizzo di librerie di sciami elettromagnetici, 
prodotti dall'interazione dei gamma con la materia. Tali librerie di sciami, prodotte con 
programmi di simulazione del CERN NASIM-GEANT, vengono impiegate da NMC per 
simulare l'energia depositata celle per cella sul calorimetro elettromagnetico. 
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4.2 Simulazione fisica 



4.2.1 I fasci K s e K L 

I K s e K L vengono generati utilizzando lo spettro di produzione 

\A 11 H P 2 n on -(8.5-E-+3.0*£) , , -(13.0-2-+3.5*£)i 

4.1 = 1.30e v po e2 ; +4.35e "o ^ 

dp dfi 4p 

dove 

[4.2] 9 S = 4.2 mrad relativo a if 5 

[4.3] ^5 = 2.4 mrad relativo a K L 

S~1 T 7" 

[4.4] p = 450 

c 

e una distribuzione di decadimento 
[4.5] d(z) = 

A 

dove 

[4.6] A = 7cr 

Lo spettro di decadimento e' ottenuto pesando lo spettro di produzione dei K con la 
distribuzione di decadimento e integrando sul volume fiduciale nel quale i decadimenti 
sono generati. 

Le particelle Ks e vengono generate entro un intervallo di energia ed entro un in- 
tervallo di vertice di decadimento tipici dell'esperimento: 



Energia min. 
Energia max. 


60 GeV 
180 GeV 


Vertice min. 
Vertice max. 


— 100 cm 
3500 cm 



Tab. 4. 2.1: Limiti valori generati 



II vertice di decadimento e' definito come distanza del punto di decadimento del K dalla 
targhetta K s ■ 

Circa 6 m oltre la targhetta Ks e' collocata l'ultima serie di collimatori di diametro 0.3 cm 
e 3.0 cm rispettivamente per Ks e Kl ■ E' percio' possibile che i gamma prodotti da un 
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K decaduto poco prima del collimatore finale riescano comunque a giungere sul rivela- 
tore attraverso il foro di apertura del collimatore stesso, e vengano quindi rivelati sul 
calorimetro elettromagnetico. Per quanto riguarda i K$ un taglio netto sul limite infe- 
riore del vertice di decadimento e' effettuato, oltre che dall'ultima serie di collimatori, 
dall'anticoincidenza K s impiegata per porre un veto sull'accettazione di quegli eventi con 
fotoni esterni all'apparato di rivelazione. 

I decadimenti di Kl distanti oltre 38 m dalla targhetta Ks non sono presi in considerazio- 
ne, in quanto l'angolo di apertura dei gamma generati da tali decadimenti non consente 
ai gamma di giungere oltre il foro del calorimetro elettromagnetico. 

II fascio di Ks , distante sul piano perpendicolare alia beam-line di 7.2 cm, all'altezza 
della targhetta Ks , dal fascio Kl , viene deflesso, attraverso il collimatore Ks , verso i 
rivelatori con un angolo di 0.6 mrad dal fascio Kl . La divergenza dei fasci Ks e Kl , 
seppur minima, e' tenuta in considerazione dal programma di simulazione. 



4.2.2 I modi di decadimento dei K 

II programma Montecarlo e' in grado di simulare, attraverso banche dati dedicate, 15 
diversi modi di decadimento del K , oltre che 3 differenti modi di decadimento del n . 
I modi di decadimento di K s e K L : 



particella 



decadimento 



Kc 



K s 
K s 



71 71 

7r°7r° 



7r ± e T u 

7T ± /X T Z/ 

7T 7r 7r° 

7T + 7T~ 

7r°7r° 

7t + 7t~e + e~~ 
7r°e + e~ 

7r 7T°7 



Kr 



K L 
K l 
K l 
K l 
K l 
K l 
K l 
K l 
K l 
K l 
K l 
K l 
K l 



Tab.4.2.2: Modi di decadimento NMC K s , K L 
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I modi di decadimento di tt : 



particella 


decadimento 


7T 


7r ± — > /j^v 
7T° — > 77 
7T° — > e + e~7 



Tab.4.2.3: Modi di decadimento NMC vr 

4.3 I rivelatori 

4.3.1 II sistema di etichettatura dei Ks 

Se la particella viene creata da un protone primario incidente sulla targhetta Ks , il 
sistema etichetta tale evento come K s ■ 

4.3.2 Le anticoincidenze e i contatori di veto 

Le anticoincidenze per i Ks e Kl , utilizzate per sopprimere i decadimenti con gamma al di 
fuori della traiettoria dei rivelatori, sono costituite da contatori a scintillazione preceduti 
da un convertitore metallico di Piombo. 

La simulazione di tali anticoincidenze avviene attraverso un calcolo probabilistico della 
conversione dei gamma nel metallo. Sono inoltre settabili differenti efficienze dell'apparato 
in funzione di 13 differenti range di energia. 

4.3.3 Lo spettrometro magnetico 

Quattro camere a fill costituiscono, insieme al magnete centrale, lo spettrometro magne- 
tico. 

La simulazione della risoluzione spaziale di tali camere e' ottenuta attraverso uno smear- 
ing gaussiano con 

[4.7] a = 103 nm 

delle coordinate trasverse effettive. 

4.3.4 L'odoscopio carico 

L'odoscopio per la rivelazione dei decadimenti carichi, necessario per una fase di pre- 
trigger, e' costituito da due semipiani verticali adiacenti di materiale a scintillazione. La 
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simulazione avviene attraverso una semplice rivelazione del punto di impatto del 7r ± sul 
corrispondente semipiano. 

4.3.5 II calorimetro adronico 

II calorimetro adronico rivela muoni /r* 1 , pioni carichi 7r =l= e protoni p. 

La simulazione di tale calorimetro fornisce l'energia del pione o del protone con un 
errore di tipo gaussiano dato da: 

dove 

[4.9] d = 0.065 

La risoluzione in posizione e' data da uno smearing gaussiano della posizione reale, 
usando una larghezza funzione dell'energia depositata: 



[4.10] a. 



x,y 



P [GeV] 



dove 

[4.11] c= 15 cm 

4.3.6 II calorimetro elettromagnetico 

II calorimetro elettromagnetico rivela gamma, elettroni e pioni, frutto del decadimento di 
K s eK L . 

La simulazione del calorimetro fornisce l'energia depositata sulle singole celle in accordo 
con la risoluzione del rivelatore: 

[4.i2] = (^) 2 + b> 

dove 

[4.13] a = 0.035 

[4.14] b = 0.005 

La risoluzione in posizione e' simulata da uno smearing gaussiano della posizione reale, 
usando una larghezza funzione dell'energia depositata: 
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[4.16] c x = 0.42 cm 



[4.17] a« = , ^ ± 0.06 cm 

[4.18] c x = 0.43 cm 

4.3.7 L'odoscopio neutro 

L'odoscopio neutro e' diviso, sia sul piano verticale, sia sul piano orizzontale, in 64 strisce 
di materiale scintillatore. La rivelazione dell'evento nell'odoscopio e' necessaria a dare lo 
start alia catena di Trigger neutro. 

4.3.8 II contatore di veto per i muoni 

II contatore di veto per i muoni consiste di tre piani di strisce di scintillatore. 

La simulazione di tale contatore tiene conto sia di un eventuale muone rivelato, sia della 
sua posizione. 



4.4 Sciami elettromagnetici 

Attraverso il programma NASIM-GEANT del CERN, vengono simulati accuratamente 
gli sciami elettromagnetici prodotti dall'interazione dei gamma con la materia. 
La probabilita' di produzione di coppie e + e _ da parte di un gamma e' data da: 

[4.19] P(x) = 1 - e-^ ax) 

dove x e' la distanza percorsa nel mezzo, pi la densita' di centri di interazione del materiale, 
a la sezione d'urto per la produzione di coppie. 

Per quanto riguarda il calorimetro elettromagnetico, attraverso il suddetto programma 
di simulazione, vengono create delle banche dati di sciami elettromagnetici sviluppati nel 
Krypton liquido, in presenza di criostato, elettrodi di rame disposti a zig-zag, separa- 
tori e distanziatori collocati tra i vari elettrodi. In corrispondenza di diverse energie e 
punti di impatto dei gamma, esistono numerosi tipi di sciami, di modo da poter simulare 
realisticamente la deposizione di energia sul calorimetro elettromagnetico. 

Per una rapida analisi dei dati prodotti con metodo Montecarlo, e'stato necessario 
ridurre la dimensione delle singole shower-box ad 11 x 11 celle, inserendo un fattore di 
correzione per tenere conto dell'energia persa con questo sistema di approssiamzione. II 
fattore di correzione stimato , ottenuto dal raffronto dei dati generati e approssimati, 
risulta essere pari a: 

[4.20] a = 1.091 
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I risultati, uniti ad una notevole velocizzazione del programma, hanno dato ragione a tale 
scelta. 

4.5 Trigger 

Sempre all'interno del programma Montecarlo NMC sono presenti le simulazioni di alcuni 
dei livelli di trigger previsti dall'esperimento NA48: 

• Trigger di livello 1 per un controllo rapido di anticoincidenze e contatori di veto 
disposti lungo la beam-line, e per una soppressione degli eventi con energia totale 
ricostruita inferiore a 20 GeV. Solo in seguito al consenso di qusto livello di trigger 
entrano in funzione i trigger di livello superiore; in questo modo si procede con una 
successiva riduzione del numero di eventi da analizzare, e con la possibilta' quindi 
di una ricostruzione piu' accurata del decadimento. 

• Trigger carico di livello 2 per la ricostruzione di energia, vertice di decadimento, 
massa invariante di 7r + 7r~ , coordinate nelle camere a fili per i decadimenti carichi. 

Presente all'interno di NMC, ma in versione alquanto semplificata, e' il Trigger neu- 
tro di livello 2, necessario alia ricostruzione di energia, baricentro dell'energia, vertice di 
decadimento e numero di gamma incidenti sul calorimetro elettromagnetico per i decadi- 
menti neutri. 

Totalmente assente e' invece il Trigger di livello 2B, [18] necessario alia distinzione, 
tra gli eventi a 5 e 6 gamma, dei decadmenti in 7r°7r 7r° da quelli in 7r°7r° con uno o due 
gamma accidentali. 
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Capitolo 5 

Decadimenti neutri di Kg e Kl 



Una accurata analisi di decadimenti K s — > 7r°7r° , K L — > tt°-k° , K L — > 7r 7r 7r prodotti con 
metodo Montecarlo ha consentito di ottenere informazioni dettagliate sulle varie quantita' 
misurabili per i tre modi di decadimento. 

Tali dati, cui e' stato necessario apportare alcune modifiche, sono stati poi elaborati 
dal simulatore del Trigger neutro, al fine di poter verificare, attraverso il raffronto delle 
quantita' ricostruite con quelle generate, l'efficienza dell'apparato stesso. 

5.1 Analisi relativist ica del decadimento 

Dalla sola conoscenza dell'energia rilasciata dai gamma sulle singole celle del calorimetro 
elettromagnetico, e' possibile ricostruire dati utili al riconoscimento dell'evento, quali 
massa invarante delle particelle in gioco, traiettorie e distribuzioni di energia dei pioni e 
dei gamma originati dal decadimento di un K . 

5.1.1 Massa invar iante 

Dalla ricostruzione della massa invariante dei n° e del K , e' ad esempio possibile distin- 
gure tra decadimenti K s — > rr°TT° , K L — > -k°-r° , e decadimenti invece K L — > 7r°7r°7r con 
soli 4 gamma sul calorimetro: in quest'ultimo caso, quand'anche fosse possibile ricostruire 
esattamente la massa invariante di due dei tre ir° , sarebbe invece impossibile ricostruire 
la massa invariante del K originario. Alio stesso modo, dalla ricostruzione della massa in- 
variante dei 7T° e del K , e' possibile distinguere tra un decadimento di tipo K s — > 7r°7r° o 
Kl — > 7r°7T con 4 gamma sul calorimetro, e un decadimento neutro di qualunque tipo con 
soli tre gamma sul calorimetro con la sovrapposizione temporale di un gamma accidentale 
sincrono con gli altri tre. 
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Per ricostruire la massa invariante dei n° e del K dall'energia depositata cella per cella, 
occorre innanzitutto riuscire a distinguere il punto d'impatto sul calorimetro e l'energia dei 
singoli gamma. Per questo occorre ricostruire esattamente le shower-box di ogni gamma. 
Tale ricostruzione risulta facile nel caso in cui le varie shower-box non siano spazialmente 
sovrapposte, e viene eseguita attraverso la ricerca del massimo di energia depositata sulle 
singole celle e il calcolo della derivata dell'energia rispetto alia posizione, al fine di trovare 
i confini della singola shower-box. La distribuzione dell'energia in funzione della posizione 
trasversa alia beam-line dello sciame elettromagnetico e' infatti nota sperimentalmente, e 
risulta dalla somma di tre esponenziali decrescenti, rappresentanti la testa, il corpo, e la 
coda dello sciame, con simmetria circolare. 

Ponendo un limite inferiore alia derivata di tale funzione, e' possibile individuare i confini 
della singola shower-box senza il rischio di conteggiare anche celle in cui la deposizione di 
energia sia puramente frutto di un basso noise interno o esterno all'apparato di rivelazione. 

Una volta ricostruite con precisione energia e punto di impatto dei gamma, si precede 
con la ricerca delle coppie di gamma che, in quanto a massa invariante, possono provenire 
dal decadimento dello stesso 7r° . Definendo un intervallo di accettazione per la massa 
invariante del 7r° in funzione dell'eflicienza e della precisione del rivelatore e del livello 
a cui si ricostruisce la quantita' fisica in esame, la massa del primo 7r° e la massa del 
secondo 7r° , nel caso di decadimenti Ks — > tt ti e Kl — > 7r°7r° , risultano essere correlate 
negativamente, con un diagramma di dispersione degli eventi di forma ellissoidale, con 
asse maggiore, inclinato di —45° rispetto alia verticale, direttamente proporzionale ai 
limiti imposti sulla massa del singolo ir° , e asse minore direttamente proporzionale alia 
risoluzione sulla misura dell'energia da parte dell'apparato di rivelazione. 
Per una buona ricostruzione degli eventi, i limiti inferiore e superiore alia massa del 
7T° (0.135 GeV) sono fissati in 0.134 GeV e 0.136 GeV rispettivamente. Da una analisi 
dei dati ottenuti, il diagramma di dispersione, per i soli eventi a 4 gamma sul calorimetro, 
risulta avere una larghezza di 383 KeV, corrispondente a 3.8 • 10~ 3 m(7r°), dello stesso 
ordine di grandezza della risoluzione in energia del calorimetro eletromagnetico 

rK11 °e 0.035 

[5 ' 1] -W = W 



5.1.2 Distribuzione di gamma e pioni 

In seguito alia identificazione delle singole particelle (pioni e gamma) prodotte dal decadi- 
mento del K , si riesce a ricostruire la traiettoria di volo virtuale dei n° (virtuale perche' 
il 7T° decade subito, in conseguenza di una vita media estremamante breve) e l'angolo di 
apertura delle coppie di gamma provenienti dal singolo 7r° . La traiettoria di volo del 
7T° viene costruita come congiungente del vertice di decadimento del K con il baricentro 
dell'energia dei due gamma provenienti dal 7r° . 
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Per quanto riguarda l'angolo di apertura dei gamma, questo e' stato calcolato nel riferi- 
mento del laboratorio, e ricostruito nel riferimento della particella tt° . 

Nel riferimento del laboratorio sono state ottenute le seguenti distribuzioni, rispettiva- 
mente per l'angolo di apertura di ciascuna coppia di gamma, e per l'angolo rispetto alia 
beam-line della traiettoria del ir° : 

Fig.5.1.1: K s — > 7r°7r° Angolo di apertura di gamma e 7r° 



Fig.5.1.2: Kl — > 7r°7r° Angolo di apertura di gamma e 7T 1 



Fig.5.1.3: K L — > 7r 7r°7r° Angolo di apertura di gamma e n° 

Per il decadimento Kl — > 7r°7r 7r° a tre corpi, l'angolo medio di apertura dei gamma 
risulta maggiore del corrispondente angolo per il decadimento Kl — > 7r°7r° a due corpi; 
i gamma, 6 invece di 4, sono mediamente meno energetici, e per la conservazione del 
quadriimpulso totale del sistema, l'angolo 9^ di apertura dei due gamma deve essere 
mediamente maggiore. La conservazione del quadriimpulso totale implica 

[5.2] m^o 2 = E^E^l - cos %) 

Sempre per il decadimento Kl — > 7r 7r°7r° , l'angolo tra beam-line e traiettoria virtuale 
del 7T° risulta minore del corrispondente angolo per il decadimento K L — > 7r°7r° ; per la 
conservazione del quadriimpulso totale del sistema (n e' il numero di rc° ): 

[5.3] m K 2 = E n oE n o(l — cos ipij) + nm^o 2 

dove ifij e' mediamente doppio dell'angolo in questione, e tenendo conto che i 37r° , seppur 
singolarmente meno energetici, sono in numero maggiore rispetto ai 7r° di un decadimento 
K L — > 7T°7r , l'angolo <pij deve essre mediamente inferiore. 

Nel riferimento della particella 7r° , l'angolo totale di apertura di ciascuna coppia di 
gamma, ottenuto attraverso una trasformazione di Lorentz delle coordinate e degli impulsi 
del laboratorio, risulta pari a 180° ± 2.7°, in accordo con un decadimento a due corpi. 
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Punto di impatto ed energia dei singoli gamma sul calorimetro risultano cosi' distribuiti: 

Fig.5.1.4: K s — > 7r°7r° Punto di impatto e Energia gamma 

I grafici si riferiscono a 10000 eventi per ciascuno dei tre modi di decadimento; per 
quanto riguarda i decadimenti Ks — > 7r°7r° occorre tenere presente che solo nel 62.4% dei 
casi almeno un gamma riesce a raggiungere il calorimetro elettromagnetico, a causa della 
presenza di una anticoincidenza anti-i^s nell'apparato. 

Fig. 5. 1.5: K L — > 7r°7r° Punto di impatto e Energia gamma 

Fig.5.1.6: K L — > 7r 7r°7r° Punto di impatto e Energia gamma 

Da una analisi dell'energia dei singoli gamma in funzione del punto d'impatto sul 
calorimetro, i gamma piu' energetici risultano fortemente concentrati verso il centro del 
calorimetro, in accordo con la conservazione del quadriimpulso totale del sistema {E 1 
maggiore implica infatti # 7 minore). La distribuzione dell'energia del gamma in funzione 
della distanza dal centro del calorimetro del punto d'impatto risulta essere: 

Fig. 5. 1.7: Ks — ► 7r°7r° Energia gamma vs. distanza dal centro 
Fig. 5. 1.8: K L — > 7r°7r° Energia gamma vs. distanza dal centro 
Fig. 5. 1.9: K L — > 7r 7r°7r° Energia gamma vs. distanza dal centro 



5.2 Analisi off-line 

L'analisi dell'evento effettuata a partire dalla conoscenza dell'energia depositata su cias- 
cuna cella del calorimetro elettromagnetico potra' essere in seguito utilizzata per una 
analisi off-line degli eventi registrati su supporto magnetico, attraverso un sistema di 
compattazione dei dati e di zero suppression [19]. 

Tale analisi consentira' di eliminare quasi totalmente gli eventi di fondo per entrambi i 
modi di decadimento, carico e neutro, e di ricostruire gli eventi utili alia misura del doppio 
rapporto R con la precisione richiesta. 



Capitolo 6. Simulazione del Trigger neutro 



36 



Capitolo 6 

Simulazione del Trigger neutro 



6.1 II rumore 

Essenziale ai fini di una simulazione realistica degli eventi e' la trattazione del rumore. 
Internamente all'apparato esistono fonti di rumore identificabili essenzialmente negli am- 
plificatori, preamplificatori e shaper collegati alle singole celle del calorimetro. Questo 
viene simulato attraverso un noise con distribuzione gaussiana a media /i — e sigma 
a = 30 MeV , corrispondente al noise massimo sperimentalmente misurato su ogni sin- 
gola cella del calorimetro; la media /i — e' frutto di un noise fluttuante nel tempo a 
media nulla. 

I noises relativi ai 128 x 128 canali del calorimetro vengono variati ad ogni ciclo di clock 
rigenerando ogni volta 128 • 128 numeri distribuiti in modo gaussiano. 

6.2 Amplificatori 

II sistema di lettura del calorimetro elettromagnetico e' costituito da due sezioni: la prima 
analogica, la seconda digitale. 

Un amplificatore amplifica il segnale proveniente da ciascuna cella del calorimetro; il 
segnale viene poi inviato ad uno shaper per formare l'impulso da analizzare. 

II segnale di ciascuna cella, formato dalla raccolta di corrente sugli elettrodi che fungono 
da condensatore, ed amplificato dall' amplificatore, e' di forma triangolare, ed e' caratte- 
rizzato da una ampiezza di 2AfiA/GeV , da un tempo di salita di 3 fis, e da un tempo di 
discesa, corrispondente al tempo di scarica del condensatore, di circa 50 fis [10]. 

La simulazione degli amplificatori richiede una banca dati contenente i fattori di ampli- 
ficazione, con massimo normalizzato ad uno, relativi ad ogni singolo canale, con la pos- 
sibilita' di tenere conto di eventuali malfunzionamenti, parziali o totali, di alcuni canali. 
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Non essendo ancora disponibile, in quanto non ancora effettuata, una mappatura com- 
pleta dei singoli canali, viene adottata una distribuzione di tipo gaussiano dei 128 x 128 
fattori di amplificazione, con media fj, — 1 e sigma a = 0.1. 

Fig. 6. 2.1: Amplificazione 

Tale a e' in accordo con l'errore introdotto dai condensatori utilizzati negli amplificatori. 

Ogni canale risulta poi essere sfasato, sempre a causa delle piccole diversita' dei vari 
componenti elettronici, rispetto al canale adiacente, con una distribuzione dei ritardi 
canale-canale di tipo gaussiano, con media /i = 0e sigma a = 1 ns, corrispondente alio 
sfasamento medio dei canali l'un l'altro. 

Fig. 6. 2. 2: Ritardi canale-canale 



6.3 Shaper 

Funzione dello shaper e' quella di sagomare il segnale per una successiva analisi. In cor- 
rispondenza ad una energia depositata sulla singola cella di 10 GeV, il massimo dell'uscita 
dello shaper e' di circa 2 V [13]. La forma d'onda del segnale in uscita dallo shaper, calco- 
lata per via teorica attraverso la funzione di trasferimento in/ out dell'apparato, e ottenuta 
da fit di dati sperimentali, e' caratterizzata da una altezza proporzionale alia derivata del 
segnale proveniente dall'amplificatore, da una larghezza a meta' altezza pari a 60 ns, e 
da un undershoot del 3% di durata 3 us [10]. 

Fig.6.3.1: Shaper: Amp. vs. Time (ns) 

Essendo lo shaper sensibile alia variazione del segnale in ingresso (rivela infatti la 
derivata di tale segnale), risponde immediatamente all'impulso amplificato proveniente 
dalla cella del calorimetro, e in modo del tutto indipendente dalla durata del segnale in 
ingresso. 

Pertanto, in caso di sovrapposizione temporale di piu' eventi sulla singola cella, lo shaper 
ad essa collegato riesce a distinguerli entro l'intervallo di durata del segnale di shaper, 
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ovvero 3 /is. In realta', considerando il fatto che l'undershoot di tale segnale e' relativa- 
mente basso (3%), la sovrapposizione di un secondo segnale con l'undershoot del primo 
non incide in maniera rilevante ai fini della rivelazione dell'evento, almeno al fine di un 
conteggio dei gamma incidenti sul calorimetro. Tali sovrapposizioni avvengono comunque 
molto raramente, essendo il rate previsto di eventi sulla singola cella inferiore ad 1 KHz. 

La sezione di shaper richiede particolare cura e attenzione nella simulazione. La strut- 
tura delle 128 x 128 celle del calorimetro, prodotta sotto forma di matrice, viene fornita 
di una terza dimensione, quella temporale. I metodi utilizzati sono diversi a seconda del 
tipo di analisi da effettuare. 

Un primo sistema, estremamente veloce ma leggermente riduttivo, consiste in una anali- 
si temporalmente limitata a 200 ns intorno al picco in energia rivelato, e nella soppressione 
dell'undershoot, con una conseguente concentrazione degli eventi, compatibilmente con i 
vari stadi di pipeline dell'elettronica, ed una riduzione dei tempi necessari all'analisi. 
Un secondo sistema, completo ed esatto, consiste invece nell'analisi del segnale dello shaper 
per tutta la sua durata temporale, tenendo cosi' conto anche di eventuali sovrapposizioni 
di due o piu' segnali. 

Nonostante questa sezione lavori ancora in modo continuo, e' possibile simularla cam- 
pionando il segnale in funzione di una successiva lettura da parte di una sezione digitale 
regolata da clock. 

La durata del campionamento, nella simulazione dell'apparato, e' quindi funzione del tipo 
di analisi da effettuare: 8 campionamenti sono sufficienti a coprire temporalmente il picco 
del segnale, mentre ne occorrono circa 125 per coprire anche tutto il periodo corrispon- 
dente all 'undershoot, con un rallentamento notevole del programma. 

La forma d'onda del segnale di shaper, nota analiticamente, viene inizialmente campio- 
nata a step di 1 ns e registrata in una banca dati, cosi' da rendere piu' veloce il sistema 
di simulazione. II calcolo attraverso la funzione di trasferimento richiede infatti tempi di 
esecuzione quattro volte superiori. 

Ad ogni step, corrispondente ad un tempo di 25 ns, viene letto il valore dell'energia 
depositata in ogni cella, convertito in Volt, moltiplicato per il fattore di amplificazione e 
sfasato del ritardo relativo al corrispondente canale, moltiplicato per la fase dovuta alia 
lettura da parte dell'elettronica, e sporcato dal rumore relativo al canale, quindi viene 
creato il vettore contenente la campionatura del segnale come formato dallo shaper, nel 
caso di energia diversa da zero, e shiftato di una posizione temporale. 
Nel caso di sovrapposizione temporale di ue o piu' segnali sulla stessa cella, alia nuova 
campionatura viene aggiunto, step per step, il corrispondente valore precedentemente regi- 
strato. Per le celle ove non sia depositata energia il processo e' notevolmente semplificato, 
e consiste in un semplice shift dei vari piani temporali della matrice di shaper, cosi' da 
poter simulare correttamente i vari stadi di pipeline. 
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Un differente sistema di simulazione, quale ad esempio una campionatura successiva 
del segnale, avrebbe richiesto tempi di esecuzione maggiore, e non avrebbe comunque 
apportato nessuna variazione agli stadi successivi. 

6.4 II clock 

La sezione digitale dell'apparato di trigger, strutturata in pipeline, e' progettata per 
funzionare con un periodo di clock di 25 ns, pari ad una frequenza di 40 MHz. Tale 
clock, totalmente indipendente dall'arrivo dell'evento sul calorimetro, non consente una 
campionatura del segnale di shaper sincrona col tempo relativo al proprio massimo, ma 
fornisce una campionatura con fase distribuita in maniera casuale uniforme in un intervallo 
di 25 ns rispetto al tempo del picco. 

In generale non e' quindi disponibile la lettura del massimo del segnale, ma un valore 
compreso tra il 75% e il 100% del massimo. Compito di una parte dell'apparato sara' poi 
quello di ricostruire, attraverso una interpolazione di tipo parabolico, il valore del picco, 
corrispondente al valore dell'energia rilasciata sul calorimetro. 

6.5 Struttura pipelined 

II rate elevato di eventi sul calorimetro richiede un sistema di lettura veloce per l'apparato 
di rivelazione. 

La realizzazione di un trigger tradizionale, che blocca cioe' l'acquisizione di qualsiasi 
evento per tutto il tempo necessario all'analisi dell'ultimo evento acquisito, avrebbe posto 
seri limiti alia velocita' di lettura e all'emcienza dell'apparato stesso. L'utilizzo di una 
struttura pipelined consente di ovviare a questo problema. 

La pipeline consiste nella suddivisione di una singola rete combinatoria in piu' blocchi, 
attraverso l'interposizione di registri [16]; ciascun blocco realizza una parte dell'operazione 
cui l'intera rete e' preposta, permettendo un'acquisizione piu' veloce dei dati, e quindi una 
minimizzazione dei tempi morti. 

E' possibile osservare che, dopo un periodo iniziale pari a n periodi di clock, se n e' 
il numero di stadi di pipeline, i risultati dell'operazione compiuta dal circuito vengono 
forniti con frequenza pari a quella del clock. 

La simulazione dei vari stadi di pipeline, ridotti al numero minimo necessario per una 
analisi dell'evento da parte del programma, avviene attraverso uno shift temporale delle 
varie quantita' memorizzate, ovvero attraverso uno shift degli indici di ciascuna variabile. 
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6.6 I sommatori analogici 

I segnali provenienti dagli amplificatori e dagli shaper connessi ad ogni singola cella del 
calorimetro elettromagnetico vengono sommati a gruppi di 8 x 2 in verticale e in orizzon- 
tale, ottenendo cosi' una doppia mappatura (nella vista x e nella vista y) concentrata del 
calorimetro. 



6.7 I digitalizzatori 

II segnale ottenuto dai sommatori analogici viene passato alia sezione digitale del Trigger 
neutro. Innanzitutto il segnale viene digitalizzato attraverso Fast ADC a 10 bit che lavo- 
rano ad una frequenza di 40 MHz. In uscita agli ADC il valore corrispondente all'energia 
rilasciata sul singolo blocchetto 8x2 celle si presenta come numero di step di digitaliz- 
zazione, piu' un livello di zero settato a 100 step di ADC. 
L'algoritmo impiegato e': 

, _ , (max — min) 

[6.1] step - 



(2 nbit - lev Zero- 1) 

[6.2] Digit = A ^o9-mod(Analog,step) + 

step 

dove max corrisponde al livello massimo di saturazione del singolo Fast ADC (100 GeV), 
min e' fissato a 0, step corrisponde all'ampiezza del singolo passo di digitalizzazione, nbit 
e' il numero di bit dei Fast ADC, levZero e' il livello di piedistallo. 

II piedistallo di 100 step di digitalizzazione e' introdotto negli ADC per poter tenere 
conto dell'undershooot dei segnali in alcuni tipi di analisi; tale livello viene invece sop- 
presso attraverso un semplice sottrattore di piedistallo (baseline restoring) 

[6.3] Digit = Digit — levZero 

[6.4] if (Digit. It. lev Zero) Digit = 

per ripristinare lo zero del segnale per la successiva analisi da parte del Trigger neutro. 

Le variabili sono trattate come Real per sfruttare al meglio le potenzialita' dei compila- 
tori Fortran, quindi velocizzare il processo, senza perdere in precisione: l'unica differenza 
infatti nella registrazione in memoria di un numero Real rispetto ad un Integer sta nel fatto 
che il primo e' caratterizzato da virgola mobile, conservando pero' un numero maggiore 
di cifre significative nella propria mantissa. 

II procedimento di digitalizzazione, come pure la somma analogica in blocchetti 8x2, 
avviene separatamente per le due viste x e y. 
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Considerando un segnale in ingresso distribuito in modo casuale uniforme tra e 
100 GeV, ed un numero di bit dell'ADC pari a 10, il rapporto tra numero digitalizzato e 
analogico 

\a <\ Digit 

Analog 

risulta avere una ditribuzione siffatta: 

Fig. 6. 7.1: Quantizzazione ADC 



con un evidente effetto di quantizzazione, presente in maniera predominante sui valori 
inferiori. 

E' lecito quindi aspettarsi una perdita in precisione sulle basse energie, poco rilevante 
comunque per Tidentincazione dell'evento e per la successiva analisi effettuata dal trigger. 



6.8 I filtri 

II segnale digitalizzato viene inviato ad un sistema di analisi, costituito dal filtro tempo- 
rale e dalla maschera spaziale. 

Entrambi questi sistemi lavorano, indipendentemente sulle due proiezioni, sulle somme 
digitali dei singoli canali (d'ora in poi sara' utilizzata la notazione canale anziche' bloc- 
chetto di 8 x 2 celle). 

II filtro temporale abilita il singolo canale, quello temporalmente successivo e i due 
precedenti, nel caso in cui il valore digitalizzato risulti superiore ad una soglia settata 
intorno a 0.5 GeV, ovvero poco al di sopra del valore massimo di noise raggingibile sul 
singolo canale, e corrispondente a circa 30 MeV ■ 16 = 0.48 GeV. Questo comporta 
una riduzione del rumore, abilitando soltanto i gruppi di celle che possono effettivamente 
contenere un picco in energia. 

Una ulteriore riduzione del rumore viene operata attraverso una selezione dei canali 
che devono contribuire alia successiva somma in righe e colonne del calorimetro elettro- 
magnetico. Questa selezione viene eseguita dalla maschera spaziale, che, confrontando il 
segnale proveniente da ciascun canale con una soglia prestabilita di 0.5 GeV, abilita alia 
somma esclusivamente quei canali per cui il segnale risulta essere sopra soglia, e un certo 
numero di canali adiacenti, a seconda della configurazione selezionata. 
L'utilizzo della maschera spaziale nell'esperimento NA48 e' tuttora incerto. 

La simulazione di entrambi questi apparati, non complessa dal punto di vista con- 
cettuale, richiede particolare attenzione nella ricostruzione dei vari stadi di pipeline, con 
un continuo aggiornamento dei canali abilitati alle somme in righe e colonne. 
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6.9 I sommatori in righe e colonne 

I canali abilitati dal filtro temporale e dalla maschera spaziale vengono sommati in righe 
e colonne separatamente nelle due viste x e y, ottenendo in questo modo 64 valori cor- 
rispondenti alle 64 proiezioni per ciascuna vista. 

I sommatori digitali a 12 bit eseguono ad ogni step temporale di 25 ns la somma dei 
vari canali, tenendo conto delle rispettive abilitazioni per i vari stadi di pipeline; quattro 
sono infatti le informazioni temporali necessarie al Peak-Finder per la ricostruzione di un 
picco. 

La variazione dell'abilitazione fornita dai filtri in slot temporali successive e precedenti 
quella in esame, obbliga ad eseguire la somma ogni volta per vari stadi di pipeline. 

Ricostruendo l'evoluzione temporale di una delle 64 colonne, tenuto conto di sfasamento 
e fattore di amplificazione delle singole celle, fase del clock rispetto al tempo del picco, 
filtro temporale, e raffrontando il segnale cosi' ottenuto col segnale di uno degli shaper, si 
ottiene: 

Fig.6.9.1: Colonna digitalizzata e Shaper 

L'informazione originaria sull'energia depositata su ogni singola cella e' adesso compat- 
tata in informazioni sull'energia depositata su una colonna o riga di 128 x 2 celle, con 
una inevitabile perdita in risoluzione, ma con una conseguente notevole velocizzazione del 
sistema di trigger, cosi' da rendere possibile una veloce analisi dei dati, compatibile con 
il rate di eventi sul calorimetro elettromagnetico. 

Una ricostruzione precisa dell'evento, con l'acquisizione dei dati relativi alle singole 
celle, viene effettuata alia fine della catena di trigger, nel caso in cui quest'ultimo riveli 
l'evento come candidato ad essere un decadimento K s — > 7r°7r° , K L — > 7r°7r° , e dia quindi 
l'abilitazione alia registrazione dell'evento per una successiva analisi off-line. 

6.10 II Peak-Sum 

Elemento di fondamentale importanza nell'apparato di Trigger neutro e' il Peak-Sum [16], 
progettato per la ricostruzione, separatamente nelle due viste, di energia e momenti di 
ordine 1 e 2 dell'energia rispetto alia posizione proiettata. 

II Peak-Sum elabora gli ingressi costituiti dalle 64 somme in righe e colonne attraverso 
l'impiego di singoli chips, progettati per poter lavorare in differenti modi, ciascuno dei 
quali vede quattro canali. Sedici chips sono quindi necessari per l'analisi di 64 canali. 

Gli algoritmi impiegati sono: per la ricostruzione dell'energia (momento di ordine 0): 



[6.6] 



M 0X = Y,E itl 
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per la ricostruzione dei momenti di ordine 1: 



[6.8] 



M 1X = 



E EjXj 



[6.9] 



M 1Y = 



per la ricostruzione dei momenti di ordine 2: 



[6.10] 



M 2X 



E EjXf 
EEi 



[6.11] 



M 2Y 



E EjY? 
EEi 



dove Xi e Yj sono rispettivamente il numero di colonna e riga nelle proiezioni x e y, ed Ei 
l'energia ricostruita sul rispettivo canale. 

II momento di ordine corrisponde esattamente all'energia depositata e ricostruita sulla 
singola vista. 

II momento di ordine 1 corrisponde al baricentro dell'energia depositata sul calorimetro 
elettromagnetico. 

II momento di ordine 2 viene utilizzato, unitamente ai precedenti, per il calcolo del vertice 
di decadimento della particella, attraverso la formula 



dove Z v e' la distanza del vertice di decadimento dalla targhetta K s , D la distanza di 
tale targhetta dal calorimetro elettromagnetico (122.35 m), E l'energia totale ricostruita, 
mk la massa del K (0.497 GeV). La formula per il calcolo del vertice di decadimento Zy 
e' data dalla conservazione del quadriimpulso totale del sistema, nell'approssimazione di 
piccoli angoli , e dalla conoscenza della massa invariante del K . 

Energia, baricentro dell'energia, vertice di decadimento della particella sono, insieme 
all'informazione sul numero di picchi rivelati sul calorimetro elettromagnetico, le quantita' 
fisiche dall'analisi delle quali e' possibile ricostruire il tipo di decadimento, e rimuovere 
quindi il fondo dovuto ad eventi accidentali e ai decadimenti Kl — > 7r 7r°7r° . 

Le somma per i momenti di ordine 0,1,2 dell'energia rispetto alia posizione proiettata 
vengono effettuate da sommatori floating point [15] a 9+1 bit di mantissa, ed esponente di 
3,4,5 bit rispettivamente per i tre momenti, per garantire una non saturazione del risultato. 
La registrazione dei soli 9 bit successivi al primo 1 nel numero in formato binario consente 



[6.12] 




Capitolo 6. Simulazione del Trigger neutro 



44 



di recuperare una cifra significativa del numero binario in questione, lasciando pressoche' 
inalterati i circuiti sommatori. 

II procedimento con cui vengono effettuate le somme all'interno del Peak-Sum si avvale 
di chip appositamente adibiti alia ricombinazione dei risultati provenienti dalle otto schede 
per ciascuna vista che compongono l'apparato, ciascuna delle quali comprensiva di due 
chips per il trattamento di quattro canali ciascuno. 

La simulazione di tale struttura, lavorando i computer su numeri a 32 o 64 bit, richiede 
una elaborata conversione del numero digitale a 12 bit in un numero floating point con 
mantissa di 10 bit, ed esponente di 3,4,5 bit a seconda dell'ordine del momento ricostruito. 
Per questa conversione sono create due apposite procedure, impiegabili a seconda del tipo 
di calcolatore e di compilatore utilizzati. 

Un primo e piu' semplice sistema consiste nel calcolo del numero di bit necessari in 
formato binario a contenere il numero in esame, attraverso il successivo confronto col valore 
2™ — 1, corrispondente al massimo numero registrabile in n bit. 1 numero digitale viene 
quindi arrotondato con un numero di cifre significative pari a 10 attraverso l'algoritmo: 



nel caso in cui n sia maggiore di 10, e lasciato inalterato nel caso in cui n sia minore di 
10. I numeri cosi' trattati vengono nuovamente registrati in memoria come Real. 

II secondo metodo per la conversione di un numero qualsiasi in floating point consiste 
nella costruzione di apposite maschere, di 26 o 58 bit a seconda del calcolatore e del com- 
pilatore utilizzati, dove un solo bit in posizione n e' uguale a 1, ed i restanti sono nulli. 
II numero digitale, registrato sotto forma di Real a 32 (64) bit, viene confrontato con le 
26 (58) maschere per identificare, attraverso operazioni logiche di tipo AND e OR, la 
posizione del primo bit diverso da zero nella mantissa, e quindi ricostruire la mantissa del 
numero in esame. L'esponente e' invece fornito direttamente dal calcolo del logaritmo in 
base 2 del numero. 

I calcolatori lavorano internamente in base 16, e l'incremento (diminuzione) di una unita' 
dell 'esponente corrisponde quindi ad uno shift verso sinistra (destra) della virgola di quat- 
tro posizioni; non potendo eccedere il numero di bit previsti, ne' perdere in cifre significa- 
tive, il compilatore procede ad una riallocazione dei bit (normalizzazione), in modo tale 
da riportare il primo 1 in una delle prime quattro posizioni, e all'mcremento o diminuzione 
di una unita' dell'esponente. Considerando che i compilatori Fortran registrano i numeri 
in formato Real in modo tale da disporre il primo 1 in uno dei primi quattro bit, sono 
sufficienti solo quattro raffronti per ricostruire esattamente la mantissa. 
Questo procedimento consente di estrarre dal numero Real nella memoria del calcolatore 
i primi 9 bit successivi al primo bit diverso da zero. 

La conversione da Real (26 bit di mantissa + 6 bit di esponente) a floating point (10 
bit di mantissa + 3,4,5 bit di esponente) viene effettuata dopo ogni somma in uscita da 



[6.13] 
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ciascuno dei 16 chips, indipendentemente per i momenti dell'energia di ordine 0,1,2. 



Lo stesso chip [16]che provvede alia realizzazione del Peak-Sum puo' essere programmato 
per la ricerca dei picchi di energia rilasciata sul calorimetro elettromagnetico, al fine di 
ricostruire il numero di gamma incidenti sul rivelatore. 

Perche' un picco venga rivelato come tale occorre che questo risulti dalle analisi com- 
binate spaziale e temporale. 

L'algoritmo di ricerca spaziale sui 64 canali rivela un picco in Xj(0) nel caso in cui: 



II valore in basso a destra si riferisce al numero del canale (riga o colonna), mentre la 
quantita' tra parentesi rappresenta il tempo dell'evento campionato ogni 25 ns. 
L'algoritmo di ricerca temporale sui 64 canali rivela un picco in Xj(0) se: 



Altra condizione imposta per l'identificazione di un picco su un canale e' che nella slot 
temporale precedente non fosse presente alcun picco nei due canali adiacenti. Questa 
condizione e' dovuta al fatto che una serie di gamma temporalmente separati di 25 ns, e 
spazialmente contigui di cella, paralizzerebbero l'apparato di trigger, non consentendo 
quindi l'analisi di nessuno di questi eventi; in un caso del genere infatti esisterebbe 
l'abilitazione dell'algoritmo di ricerca temporale, ma verrebbe a mancare con buona prob- 
abilita' l'abilitazione da parte dell'algoritmo di ricerca spaziale. 

La simulazione software del sistema Peak-Finder si avvale, per la ricerca dei picchi, dei 
dati registrati nei successivi stadi di pipeline, limitabili ad un numero minimo di cinque 
per una accurata simulazione dell'apparato. 

Le informazioni sul numero di picchi e sulla loro posizione, vengono registrate in appo- 
site banche dati ad accesso rapido. 

6.12 La temporizzazione dell'evento 

Lo stesso chip impiegato per la ricostruzione dei momenti di ordine 0,1,2 dell'energia 
rispetto alia posizione proiettata e per la ricerca dei picchi e' programmabile per la ri- 
costruzione del tempo del picco in sottoslot di 3.125 ns, intervallo corrispondente a 1/8 
del periodo di clock. 



6.11 II Peak-Finder 



[6.14] 
[6.15] 



Xi(0) >X i+1 (0) 
Xi(0) >X,_ 1 (0) 



[6.16] 
[6.17] 



Xi(0) > Xi(+1) 
Xi(0) > Xi(-l) > Xi{-2) 
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Attraverso una interpolazione rettilinea dei valori campionati ogni 25 ns, l'analisi pro- 
cede con una ricerca della sottoslot temporale in cui cade il tempo di riferimento. Tale 
tempo di riferimento corrisponde, a seconda della scelta programmabile via VME, al 
tempo in cui il segnale raggiunge rispettivamente |, \ + ^, \ + |, \ + | + dell'altezza 
del massimo campionato, per una migliore precisione nella temporizzazione del segnale. 

In funzione della differenza tra i tempi relativi ai vari picchi rivelati, il Trigger Supervisor 
dara' poi l'abilitazione all'acquisizione dell'evento. 

Da una analisi effettuata su 30000 eventi con gamma sincroni sul calorimetro elettro- 
magnetico, emerge la seguente distribuzione di distanza temporale massima tra i tempi 
relativi ai vari picchi, espressa in numero di sottoslot da 3.125 ns ciascuna: 

Fig. 6. 12.1: Distribuzione A T dei gamma 

Tali risultati suggeriscono quindi un intervallo temporale di accettazione per eventi 
buoni corrispondente a (0, 3) sottoslot di differenza massima tra i tempi relativi ai vari 
gamma. 

6.13 Look-up table 

I moduli Look-up table sono costituiti essenzialmente da memorie programmabili, e per- 
mettono l'elaborazione dei risultati forniti dal Peak-Sum e dal Peak-Finder per il calcolo 
delle quantita' richieste dal sistema di trigger [17]. 

Le funzioni realizzate da questi moduli riguardano l'eleborazione dei dati relativi alle 
proiezioni x e y, e consentono di determinare: 

• tempo esatto dell'evento sul calorimetro elettromagnetico 

• energia totale depositata sul calorimetro elettromagnetico 

• posizione del baricentro dell 'energia sul calorimetro 

• posizione del vertice di decadimento della particella 

• numero di picchi sul calorimetro elettromagnetico 

La ricostruzione del valore dell'energia depositata sul calorimetro elettromagnetico 
avviene attraverso una interpolazione parabolica del valore del massimo e dei due va- 
lori temporalmente contigui. 
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6.14 II Trigger Supervisor 

Una serie di tagli sulle quantita' calcolate viene effettuata dal Trigger Supervisor, in modo 
da ridurre il rate di eventi candidati ad essere registrati su supporto magnetico. 

I dati generati da NMC e quelli ricostruiti dal trigger relativi agli eventi analizzati 
vengono registrati su disco, per poter poi calcolare l'efficienza dell'apparato sperimentale 
in corrispondenza di diversi tagli effettuati. 
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Capitolo 7 
Test e risultati 

7.1 Campionamento del segnale ed errori di arroton- 
damento 

Attraverso la simulazione dettagliata di tutto l'apparato sperimentale, costituito da rive- 
latore piu' trigger, e' possibile individuare quali, tra le varie fonti di arrotondamento e 
imprecisione introdotte, apportino un maggiore errore nella determinazione delle quantita' 
fisiche finali. 

Fonti di approssimazione presenti all'interno dell'apparato di trigger sono: 

• il campionamento del segnale effettuato in maniera non sincrona col tempo del picco 
in energia sul calorimetro elettromagnetico 

• la differenze nei fattori di amplificazione dei singoli canali 

• lo sfasamento relativo canale-canale 

• la compattazione delle celle in blocchetti 8x2 

• la digitalizzazione dell'informazione 

• la possibilita' di saturazione dei digitalizzatori ai vari livelli 

• il sistema di filtri con livello minimo di accettazione (soglia) 

• la somma digitale in righe e colonne di 128 x 2 celle 

• l'analisi nelle due proiezioni separate anziche' bidimensionale e di tutte le celle del 
calorimetro senza accorpamenti 
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Da numerosi test effettuati sull'apparato attraverso la sua simulazione software emer- 
gono, in relazione alle diverse sezioni del Trigger neutro, i seguenti errori relativi medi 
sulla ricostruzione dell'energia totale, tutti distribuiti in maniera gaussiana, il cui effetto 
e' concorde nell'apportare una minorazione dell'energia ricostruita dal trigger. 

I risultati sono ottenuti da successive analisi, in ciascuna delle quali vengono annullate 
tutte le possibili fonti di errore eccetto una. 



fattore di amplificazione 


-0.004 


± 


0.043 


sfasamento canale-canale 


-0.002 


± 


0.005 


fase clock-picco 


-0.102 


± 


0.075 


ADC 10 bit 


-0.005 


± 


0.002 


filtro temporale 0.5 GeV 


-0.027 


± 


0.013 


saturazione 


-0.000 


± 


0.002 


floating point 


-0.001 


± 


0.003 


accorpamento celle 


-0.000 


± 


0.000 



Tab. 7.1.1: Errori ed arrotondamenti 



11 raggruppamento di piu' celle in gruppi non comporta, al fine di una esatta ri- 
costruzione dell'energia, alcun errore. 

Escludendo i moduli Look-up table, la ricostruzione dell'energia totale depositata sul 
calorimetro elettromagnetico avviene con un errore relativo distribuito secondo una gaus- 
siana con media \i = —0.132, a = 0.147, e troncata a 0. L'energia ricostruita dal Peak-Sum 
risulta quindi inferiore all'energia effettiva mediamente del 13%. 

Solo in seguito ai moduli Look-up table, attraverso una ricostruzione dell'energia to- 
tale ottenuta con interpolazione parabolica, l'errore medio relativo sull'energia ricostru- 
ita rispetto all'energia effettivamente depositata sul calorimetro risulta pari a —0.012 ± 
0.047, con una conseguente diminuzione dell'errore di un fattore 10 rispetto al precedente 
-0.132 ±0.093. 

Le distribuzioni dell'energia e dell'errore relativo sull'energia sono riportate, separata- 
mente per i tre modi di decadimento Kg — > 7r°7r° , Kl — > 7r°7r° , Kl — > 7r°7r 7r° , nelle 
prossime sezioni, insieme alle distribuzioni di numero di gamma sul calorimetro, baricen- 
tro dell'energia, vertice di decadimento, e errori relativi delle quantita' ricostruite rispetto 
alle corrispondenti quantita' generate. 

7.2 Distribuzione temporale degli eventi 

In funzione del tipo di analisi da effettuare, e' possibile, all'interno del simulatore di Trigger 
neutro, selezionare diverse configurazioni per la distribuzione temporale degli eventi. 
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E' possibile, ad esempio, tenere conto del branching ratio dei diversi modi di decadi- 
mento neutri di K$ e Kl , cosa peraltro inutile ai fini di una analisi della risposta 
dell'apparato ai tre diversi modi di decadimento K$ — > 7r°7r° , Kl — > 7r°7r° , Kl — > 
7T 7r 7r° , in quanto si verifica mediamente un decadimento Kl — > 7r°7r° ogni 200 decadi- 
ment K L -> 7r 7r°7r° . 

Constatato il fatto che la simulazione dell' undershoot non comporta alcuna modifica 
apprezzabile ai risultati, essendo presente nell'apparato di trigger un sottrattore di piedi- 
stallo (baseline restoring) per ripristinare il livello di zero dopo gli ADC, gli eventi vengono 
distribuiti in maniera tale da risultare temporalmente separati in media di 10 slot tem- 
porali da 25 ns, tempo sufficiente per garantire una non sovrapposizione dei segnali di 
shaper di due eventi successivi, ed una corretta simulazione del flusso di dati attraverso i 
vari stadi di pipeline. 

La distribuzione temporale degli eventi puo' essere allargata o ristretta a piacere, at- 
traverso la modifica di alcuni parametri. 

I test di cui sono riportati i risultati sono effettuati senza sovrapposizione di eventi sulla 
stessa cella del calorimetro elettromagnetico, fatto che peraltro si verifica con probabilita' 
pressoche' nulla, dato il basso rate di fotoni sulla singola cella. 

7.3 Analisi e acquisizione dei dati 

Una prima analisi dei dati generati dal programma Montecarlo consente di ricostruire 
a livello del calorimetro elettromagnetico il numero di gamma, l'energia rilasciata sul 
calorimetro, il baricentro dell'energia, il vertice di decadimento del K . Tali valori vengono 
registrati come dati di livello 1. 

I corrispondenti valori generati, quindi non ancora elaborati e filtrati dal simulatore 
dell'apparato sperimentale, nonostante non possano essere utilizzati dal trigger ne' da 
una successiva analisi off-line perche' sconosciuti, vengono comunque registrati come dati 
di livello 0. 

In seguito i dati sono elaborati dal simulatore, e registrati in uscita dal Trigger neutro 
come dati di livello 2, e come dati di livello 3 alia fine dei moduli Look-up table. 

I grafici ed i risultati riportati in seguito terranno quindi conto della seguente notazione: 



livello 0: 


NMC 


New Montecarlo 


livello 1: 


LKR 


Liquid Krypton Calorimeter 


livello 2: 


NET 


Neutral Trigger 


livello 3: 


LUT 


Look-up Table 



Tab. 7.3.1: Notazione livelli 



Per un raffronto dei dati ricostruiti a livello 1 e 3 fra loro e con i dati generati di livello 0, 



Capitolo 7. Test e risultati 



51 



e' essenziale numerare gli eventi e poter disporre in uscita del numero dell'evento ricostru- 
ito. Per questo scopo numerosi contatori sono disposti nel simulatore in corrispondenza 
di ciascuno stadio di pipeline. 

I dati a livello 1 e 3, ovvero ricostruiti da una lettura del calorimetro elettromagnetico 
e dai moduli Look-up table rispettivamente, sono riportati senza alcun taglio effettuato 
sui valori ricostruiti, al fine di mettere in evidenza il reale funzionamento del trigger. 

Solo in seguito, nel calcolo dell'emcienza dell'apparato sperimentale, verranno effettuati 
i tagli necessari sulle varie quantita' fisiche ricostruite. 

I grafici riportati di seguito si riferiscono ad una analisi effettuata su un campione di 
10000 decadimenti per ognuno dei tre modi di decadimento K s — > 7r°7r° , K L — > tt°-k° , 
Kl — > 7r 7r 7r , al fine di poter visualizzare i risultati dei confronti tra le quantita' ai vari 
livelli. 

7.4 Energia 

Le distribuzioni in energia, per i tre modi di decadimento analizzati, sono graficate con 
l'energia totale, espressa in GeV, sulle ascisse, e il numero di eevnti con energia com- 
presa neH'imtervallo di ±0.5 GeV attorno al valore del singolo step dell'istogramma sulle 
ordinate, e si riferiscono nell'ordine a: 

• Energia generata dal programma Montecarlo (NMC): ENERG 

• Energia ricostruita da Look-up Table (LUT): ENERG 3 

I grafici di raffronto si riferiscono al plot dei valori generati o ricostruiti ai diversi livelli 
per eventi corrispondenti; il primo indice si riferisce al valore riportato sulle ascisse, il 
secondo indice al valore riportato sulle ordinate del grafico. 

Fig.7.4.1: K s — > 7r°7r° Energia generata NMC e ricostruita LUT 



Fig.7.4.2: K s — > 7r°7r° Energia ricostruita LUT vs. generata NMC 

Gli eventi che, graflcati, si trovano in prossimita' della retta a 45° sono quegli eventi 
per cui tutta l'energia del K s viene depositata sul calorimetro elettromagnetico. I punti 
invece che si trovano al di sotto della retta a 45° sono riferiti a quegli eventi per cui uno o 
piu' fotoni non giungono sul calorimetro: la deposizione di energia risulta quindi inferiore 
a quella posseduta dalla particella originaria. 
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Fig.7.4.3: K L — ► 7r°7r° Energia generata NMC e ricostruita LUT 



Fig. 7.4.4: K L — > 7r u 7r u Energia ricostruita LUT vs. generata NMC 

La differente distribuzione in energia per i decadimenti Ks — > 7r°7r° da una parte, e 
K L — > ir°ir° , kltpiz dall'altra, e' dovuta essenzialmente alia configurazione dell'apparato 
sperimentale e al sistema di produzione dei due fasci Ks e Kl rispettivamente. 

Fig.7.4.5: K L — > 7r 7r°7r° Energia generata NMC e ricostruita LUT 



Fig. 7. 4. 6: Kl — > 7r 7r°7r° Energia ricostruita LUT vs. generata NMC 

Dall'analisi della distribuzione dell'energia ricostruita a livello 3 (LUT), e dal successivo 
rafTronto con l'energia generata a livello (NMC) per i tre modi di decadimento, e' 
possibile notare come per un vasto numero di eventi l'energia ricostruita sia inferiore a 
quella generata; questo effetto e' dovuto agli eventi per cui uno o piu' gamma non giungono 
sul calorimetro elettromagnetico, senza quindi che possa essere rivelata la loro energia. 

Un successivo rafTronto tra l'energia ricostruita a livello 1 (LKR) sulle ascisse e a livello 

3 (LUT) sulle ordinate mette invece in risalto la precisione dell'apparato di Trigger e 

Look-up table nella ricostruzione dell'energa totale depositata (grafici di sinistra). 

Attraverso il calcolo 
^ Energia(3) — Energia(l) 

Energia(l) 

e' possibile stimare l'errore relativo nella ricostruzione dell'energia totale da parte del 
Trigger neutro (grafici di destra). 

Fig. 7.4. 7: K s — »■ 7r°7r° Energia 1-3 e errore relativo 
K s -> 7r°7r° errore = -0.017 ± 0.038 



Fig. 7.4. 8: Kl — ► 7r°7r° Energia 1-3 e errore relativo 
K L -> 7r°7r° errore = -0.013 ± 0.046 



Fig. 7.4. 9: Kl — ► 7r 7r°7r° Energia 1-3 e errore relativo 
K L -> 7r 7r°7r° errore = -0.012 ± 0.049 
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7.5 Numero di gamma 

Mentre i gamma generati sono ovviamente 4 per i decadimenti K$ — > 7r°7r° , Kl — > 7r°7r° , 
e 6 per i decadimenti /Tl — > 7r 7r 7r° , i gamma che in realta' colpiscono il calorimetro 
elettromagnetico possono essere in numero inferiore, a causa di quei gamma che hanno 
traiettorie esterne al rivelatore e quindi anticoincisi, o che, poco energetici, non vengono 
rivelati. 

Per ognuno dei tre modi di decadimento sono riportate, nell'ordine, le seguenti distri- 
buzioni: 

• Numero di gamma che giungono sul calorimetro elettromagnetico: GAMMA 

• Numero di gamma ricostruiti sul calorimetro elettromagnetico da una analisi cella 
per cella, con una soglia per l'individuazione del picco fissata a 0.5 GeV ed un algo- 
ritmo di ricerca che prevede il confronto di ogni singola cella con le otto circostanti: 
GAMMA 1 

• Numero di gamma ricostruiti sul calorimetro dall'analisi nella proiezione x, senza 
alcun taglio effettuato, delle colonne gia' accorpate a due a due: GAMMA 1 X, 
in modo da rendere questa distribuzione confrontabile con quella ottenuta a livello 
3. La soglia per l'individuazione di un picco in una colonna e' fissata a 0.5 GeV, 
valore corrispondente alia soglia del filtro temporale, e l'algoritmo di ricerca spaziale 
prevede il confronto di ogni colonna con le due adiacenti (per le due colonne esterne 
il confronto avviene con uno zero) 

• Differenza tra il numero di picchi rivelati dal trigger (livello 3) e il numero di gamma 
ricostruiti sul calorimetro elettromagnetico dall'analisi delle celle gia' accorpate in 
colonne (livello 1) nella proiezione x: GAMMA 3-1 X 

Fig. 7. 5.1: K s — > ir°ir° Numero gamma generati e ricostruiti cella per cella 



Fig. 7. 5. 2: Ks — ► n n Numero gamma in proiezione x e errore di ricostruzione 

Per i decadimenti Ks — > 7r°7r° , il numero di eventi in cui almeno un gamma giunge 
sul calorimetro e' pari al 62.4% del totale dei decadimenti K s — > n n° prodotti. Tale 
riduzione e' dovuta essenzialmente ai collimatori e alia anticoincidenza anti-i^ , impiegata 
per sopprimere i decadimenti Ks — > 7r°7r° con fotoni esterni al cono di accettanza del 
rivelatore, e quindi anticoincisi. 
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Fig. 7.5. 3: K L — > 7r°7r° Numero gamma generati e ricostruiti cella per cella 



Fig. 7. 5. 4: K L — ► 7r°7r° Numero gamma in proiezione x e errore di ricostruzione 



Fig. 7. 5. 5: K L — ► 7r 7r°7r° Numero gamma generati e ricostruiti cella per cella 



Fig. 7. 5. 6: Kl — ► 7r 7r°7r° Numero gamma in proiezione x e errore di ricostruzione 



Per i decadimenti Kl — > 7r°7r° e Kl — > 7r 7r 7r il numero di eventi in cui almeno un 
gamma giunge sul calorimetro elettromagnetico e' pari a circa il 99% del totale dei decadi- 
menti prodotti all'interno della zona fiduciale, in quanto non esiste un anticoincidenza 
qhH-Kl subito oltre la targhetta Kl ■ 

E' evidente come nelle analisi per colonne di 128 x 2 celle, effettuate a livello del 
calorimetro elettromagnetico e a livello del Trigger neutro, si abbia una diminuzione del 
numero di eventi rivelati a 4(6) gamma, rispetto al numero di eventi generati a 4(6) 
gamma. Tale effetto e' dovuto a eventi in cui piu' gamma colpiscono il calorimetro elet- 
tromagnetico in un intervallo spaziale orizzontale di 4 cm, corrispondente alia larghezza 
di 2 celle, ovvero di una colonna, venendo quindi rivelati come un singolo picco all'interno 
della medesima colonna nella proiezione x. 

Un numero di gamma ricostruito superiore al numero di gamma incident i sul calorimetro 
e' invece dovuto o ad effetti di sovrapposizione di due o piu' sciami elettromagnetici 
di differenti gamma, che generano quindi dei punti di sella nella distribuzione spaziale 
dell'energia, o ad effetti di bordo, consistenti nel confronto delle due colonne esterne con 
la rispettiva colonna interna adiacente da una parte, con uno zero dall'altra, algoritmo 
questo utilizzato dal Peak-Finder per l'individuazione dei picchi. 

Per la proiezione y i risultati e i grafici sono pressoche' identici. 
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7.6 Baricentro dell'energia 

II baricentro dell'energia depositata sul calorimetro elettromagnetico e' calcolato, sepa- 
ratamente per le proiezioni x e y, attraverso la formula 

[7.2] M 1X -- EElXl 



[7.3] M 



1Y 



come momento di ordine 1 dell'energia rispetto alia posizione proiettata. 

I plot del baricentro dell'energia, espresso in cm, per i tre modi di decadimento, si 
riferiscono, nell'ordine, a: 

• Plot cartesiano del baricentro, calcolato dal punto di impatto dei gamma generati 
sul calorimetro elettromagnetico: BARIC 

• Distribuzione della distanza del baricentro dell'energia, calcolato a partire dalla 
conoscenza del punto d'impatto dei singoli gamma sul calorimetro, dal centro del 
calorimetro elettromagnetico: BARIC 

• Plot cartesiano del baricentro dell'energia depositata sul calorimetro ricostruito dal 
Trigger neutro: BARIC 3 

• Distribuzione della distanza del baricentro dell'energia, ricostruito dal Trigger neu- 
tro, dal centro del calorimetro elettromagnetico: BARIC 3 

Fig. 7.6.1: K s — > 7r°7r° Baricentro generato 



Fig. 7.6. 2: K s — > 7r°7r° Baricentro ricostruito 



Nella ricostruzione del baricentro dell'energia depositata sul calorimetro elettromagne- 
tico, occorre tenere presente che la posizone cartesiana dei gamma non e' nota, bensi' se 
ne conoscono le due posizioni proiettate sulle viste x e y, con una precisione dell'ordine 
di 4 cm, corrispondenti alia larghezza di ciascun canale. 



Fig. 7. 6. 3: K L —*■ 7r°7r° Baricentro generato 
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Fig. 7.6.4: K L — > 7r°7r° Baricentro ricostruito 



Fig. 7.6. 5: Kl — > 7r 7r°7r° Baricentro generato 



Fig. 7.6. 6: K L — > 7r 7r 7T Baricentro ricostruito 



Attraverso il calcolo 
,^ ^ Baricentro(3) — Baricentro(l) 

Baricentro(l) 

e' possibile stimare l'errore relativo sul baricentro dell'energia ricostruito dal Trigger neu- 
tro in funzione dello stesso ricostruito dall'analisi delle singole celle del calorimetro elet- 
tromagnetico accorpate in righe e colonne 128 x 2. 

L'errore nella ricostruzione risulta essere, per i tre modi di decadimento: BARIC 3/1. 

Fig. 7. 6. 7: K s — > 7r°7r° Errore sulla ricostruzione del baricentro 
K s -> 7r°7r° errore = 0.004 ± 0.083 



Fig. 7. 6. 8: K L — > 7r°7r° Errore sulla ricostruzione del baricentro 
K L -> ttV errore = 0.018 ± 0. 101 



Fig. 7. 6. 9: K L — > 7r 7r°7r° Errore sulla ricostruzione del baricentro 
K L -> 3vr errore = 0.021 ± 0.196 
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7.7 Vertice di decadimento della particella 

I plot per i tre modi di decadimento del vertice di decadimento della particella K s o K L , 
espresso in cm, si riferiscono nell'ordine a: 

• Distribuzione del vertice di decadimento del K generato da NMC nella zona fidu- 
ciale: VERTX 

• Raffronto tra vertice generato da NMC a livello (sulle ascisse) e ricostruito sul 
calorimetro elettromagnetico a livello 1 (sulle ordinate): VERTX 0-1 

• Distribuzione del vertice di decadimento del K ricostruito dal Trigger neutro senza 
tagli effettuati: VERTX 3 

• Raffronto tra vertice generato da NMC a livello (sulle ascisse) e vertice ricostruito 
dal trigger a livello 3 (sulle ordinate): VERTX 0-3 

Fig. 7. 7.1: K s — ► 7r°7r° Vertice generato e raffronto NMC-LKR 



Fig.7.7.2: K s — > 7r°7r° Vertice ricostruito e raffronto NMC-NET 

II vertice generato e ricostruito per i decadimenti Ks presenta un taglio inferiore ab- 
bastanza netto in corrispondenza dell'anticoincidenza e dell'ultimo collimatore Ks , a 
6 m dalla targhetta K s ; attraverso il foro del collimatore riescono comunque a passare 
i gamma prodotti da alcuni decadimenti di Ks precedenti al collimatore stesso, il che 
giustifica una distribuzione non nulla prima dei 6 m. 



Fig.7.7.3: K L — ► 7r°7r° Vertice generato e raffronto NMC-LKR 



Fig. 7. 7.4: K L — »■ 7r°7r° Vertice ricostruito e raffronto NMC-NET 
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Fig.7.7.5: K L -> 7r 7r°7r° Vertice generato e raffronto NMC-LKR 



Fig. 7. 7.6: K L -> 7r 7r°7r° Vertice ricostruito e raffronto NMC-NET 

II vertice generato e ricostruito per i decadimenti K L — > 7r°7r° e K L — > 7r 7r 7r presenta 
una distribuzione non nulla per distanze dall'origine del sistema di riferimento inferiori a 
4.8 to, dove e' posizionato l'ultimo collimatore Kl ; attraverso il foro del collimatore ri- 
escono a passare i gamma prodotti da decadimenti di alcuni precedenti al collimatore 
stesso, il che giustifica una distribuzione non nulla prima di 4.8 m. 

Dall'analisi delle distribuzioni del vertice di decadimento ricostruito a livello 3 (NET) 
e dal successivo raffronto col vertice generato a livello (NMV) , per i tre modi di decadi- 
mento, e' possibile notare come per un vasto numero di eventi il vertice di decadimento 
ricostruito risulti superiore a quello generato. Questo e' dovuto al fatto che per gli eventi 
in cui uno o piu' gamma non giungono sul calorimetro, il vertice di decadimento, in ac- 
cordo con le previsioni, viene ricostruito ad una distanza maggiore dalla targhetta K$ , 
origine del sistema di riferimento adottato. 

La formula utilizzata infatti per il calcolo del vertice di decadimento della particella e': 

[7.5] Z v = D - —^{M 2X + M 2Y ) - (Mf x + M? Y ) 

dove Z v e' la distanza del vertice di decadimento dalla targhetta K s , D la distanza di 
tale targhetta dal calorimetro elettromagnetico (122.35 to), E l'energia totale ricostruita, 
TOfc la massa del K (0.497 GeV), M\xy e M 2 x,y rispettivamente i momenti dell'energia 
di ordine 1 e 2 rispetto a X, Y. 

La formula per il calcolo del vertice di decadimento Zy e' data dalla conservazione del 
quadriimpulso totale del sistema, nell'approssimazione di piccoli angoli, e dalla conoscenza 
della massa invariante del K. 

Un successivo raffronto tra il vertice di decadimento ricostruito a livello 3 (NET) e a 
livello 1 (LKR) mette in risalto la precisione dell'apparato di Trigger e Look-up table 
nella ricostruzione del vertice (grafici di sinistra). 

Attraverso il calcolo 

Vertice(3) — Vertice(l) 
L ' J Vertice(l) 

e' possibile stimare l'errore relativo nella ricostruzione del vertice di decadimento da parte 
dell'elettronica (grafici di destra). 
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Fig.7.7.7: K s 

K S -»• 7l°7l° 



7T 



7T Vertice 1-3 e errore relativo 
errore = 0.027 ±0.041 



Fig.7.7.8: K L -» 

K 5 -> 7T 7T 



7T 



7T Vertice 1-3 e errore relativo 
errore = 0.042 ± 0.079 



Fig.7.7.9: K L -> 

X 5 -> 7r°7r° 



7r 7r 7r° Vertice 1-3 e errore relativo 
errore = 0.051 ±0.107 



In tutti i grafici e' possibile notare una forte asimmetria tendente ad una ricostruzione 
del vertice di decadimento superiore a quello effettivo, come peraltro previsto. 

7.8 Efficienza del Trigger neutro 

In funzione di diversi parametri e di diverse configurazioni dell'apparato sperimentale si 
ottengono differenti efficienze nella rivelazione dei tre modi di decadimento K s — > 7r°7r° , 



Possibili variazioni alle soglie dei filtri temporali, un differente sistema di raggruppa- 
mento delle celle in blocchetti prima, in righe e colonne poi, e diversi tagli effettuati sui 
valori ricostruiti comportano una diversa efficienza del sistema. 

Unitamente all'efficienza di rivelazione del decadimento, e' necessario considerare anche 
il numero di eventi che costituiscono il fondo Kl — > 7r 7r°7r° rivelati dal trigger come eventi 
di decadimenti Kl — > 7r°7r° , al fine di ottimizzare l'acquisizione dati in funzione di questi 
due parametri. 

Un intervallo di accettazione troppo ampio per il vertice di decadimento ricostruito, 
mentre da una parte consentirebbe una elevata efficienza nella rivelazione dei decadimenti 
Ks — > 7r°7r° e Kl — > 7r°7r° puliti, dall'altra consentirebbe a parecchi decadimenti Kl — > 
7r°7T 7r di essere rivelati come decadimenti Kl — > 7r°7r° . 

Dal momento che il branching ratio dei decadimenti Kl — > 7r 7r 7r° e' 0.217, e il branching 
ratio dei decadimenti K L — > 7r°7r° e' 0.001, e quindi in natura avviene un decadimento 
Kl — > 7r°7r° ogni 210 decadimenti K L — > 7r 7r 7r° , una inefficienza nel riconoscimento dei 
decadimenti Kl — > 7r°7r°7r comporterebbe la registrazione su supporto magnetico di un 
elevato numero di eventi indesiderati. 



K L -> 7T°7r , K L -> 7r 7r 7r° . 
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La conoscenza dell'efficienza con una precisione di circa 1-10 4 e' necessaria per una 
misura del doppio rapporto 

[7 71 N(Kl — » 7T 7T ) N(Kl — > 7T + 7T~) ^ 6/ 

con una precisione di 2 • 10~ 4 . 

Detto infatti n il numero reale di decadimenti di un certo tipo, e no il numero misurato di 
tali decadimenti, dalla conoscenza dell'efficienza e relativa a quel modo di decadimento, 
e' possibile ricostruire il numero n\ di eventi stimati come: 

[7.8] ni = * 

e 

Per una accurata misura di R occorre quindi stimare con la precisione richiesta n x per 
i quattro modi di decadimento K s — > h + -k~ , K s — > 7r°7r° , Tf^ — > 7r + 7r~ , K L — > 7r°7r° , 
ovvero l'efficienza e, dal momento che Uq risulta dall'esperimento. 

Essendo la distribuzione degli eventi nell'unita' di tempo di tipo Poissoniano, caratte- 
rizzata da media u = n, con n numero di eventi, e a = y/Jl, ovvero da errore relativo 

[7.9] = ± 

il numero di decadimenti necessari per avere una precisione sul doppio rapporto dell'ordine 
di i0" 4 e' circa iO 8 . 

Nelle seguenti tabelle sono riportati i rapporti tra il numero di eventi in uscita (out), e 
il numero totale di eventi generati da NMC (tot) e in ingresso al Trigger (in), in funzione 
di differenti configurazioni possibili dell'apparato sperimentale di Trigger neutro. I valori 
f e si riferiscono rispettivamente alia accettazione o meno dell'evento. 

Gli eventi in ingresso (in) sono accettati se si presentano come decadimenti Ks — > 
Kl — > 7r°7T con 4 gamma sul calorimetro elettromagnetico (il dato e' fornito dal 
programma di generazione degli eventi, e, benche' sconosciuto nella realta', e' utile a 
questo livello per una esatta conoscenza dell'efficienza dell'apparato). 
Gli eventi in uscita (out) sono accettati se giudicati buoni dal Trigger neutro, ovvero se, 
una volta filtrati ed elaborati dall'apparato, i valori in uscita rientrano entro gli intervalli 
di accettazione richiesti. 

II rapporto outl/inl esprime direttamente l'efficienza nel conteggio degli eventi buoni, 
mentre il rapporto outO/inO esprime l'efficienza nell'eliminazione del fondo e dei decadi- 
menti non buoni. 

I rapporti incrociati outl/inO e outO/inl tengono conto della inefficienza dell'apparato di 
trigger nell'eliminazione del fondo e nella rivelazione di eventi buoni rispettivamente. 
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II rapporto K L 2/K L 3 indica ogni quanti eventi registrati costituenti il fondo K L — > 
7T 7r 7r° viene registrato un decadimento utile Kl — > 7r 7r° , tenuto conto del branching 
ratio dei due modi dei decadimento. Tale indicazione risulta estremamente importante 
per la mole di informazioni e il rate di eventi da registrare su disco. 

I valori di seguito riportati si riferiscono ad una analisi in assenza di eventi accidentali 
in funzione di differenti configurazioni e tagli effettuati per l'accettazione delle quantita' 
fisiche ricostruite. 

Configurazione: blocchetti indica il numero di celle componenti ciascun blocchetto e la 
loro configurazione, soglia e' la soglia del filtro temporale espressa in GeV. 

Tagli: gli intervalli di accettazione per le quattro quantita' fisiche ricostruite espri- 
mono il valore minimo e massimo rispettivamente per l'accettazione dell'evento in esame; 
gamma e' il numero di gamma nella singola proiezione, energia e' l'energia totale ri- 
costruita del K espressa in GeV , baricentro e' la distanza, espressa in cm, del baricentro 
dell'energia ricostruito dal centro del calorimetro elettromagnetico, vertice e' il vertice di 
decadimento della particella espresso in cm. 



blocchetti 


8x2 


soglia 


0.5 


energia 


(60,180) 


baricentro 


(0,10) 



Tab. 7.8.1: Configurazione e tagli 

I tagli su energia e baricentro dell'energia sono dedotti dall'analisi delle distribuzioni 
delle suddette quantita' generate e ricostruite. Per quanto riguarda invece numero di 
gamma e vertice di decadimento, in funzione dei differenti tagli applicati, si ottengono 
differenti valori di efficienza e quindi di rate di eventi da registrare su supporto magnetico. 
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gamma 






3,4 


2,5 


3,4 


vertice 






300,1500 


200,2000 


200,2000 




in 


out 








K s -> 7r°7r° 








0.997 


0.962 


0.890 







1 


0.003 


0.038 


0.110 


in0/tot=0.666 


1 





0.641 


0.244 


0.521 


inl/tot=0.334 


1 


1 


0.359 


0.756 


0.479 




tot 


1 


0.112 


0.278 


0.167 


K L 7r°7r° 








0.996 


0.972 


0.992 







1 


0.004 


0.028 


0.008 


in0/tot=0.603 


1 





0.854 


0.615 


0.774 


inl/tot=0.397 


1 


1 


0.146 


0.385 


0.226 




tot 


1 


0.060 


0.170 


0.095 


\/ , _0_0_0 
l\ 1 — >■ 7T 7T 7T 








0.991 


0.939 


0.985 







1 


0.009 


0.061 


0.015 


in0/tot=1.000 


1 





0.000 


0.000 


0.000 


inl/tot=0.000 


1 


1 


0.000 


0.000 


0.000 




tot 


1 


0.009 


0.061 


0.015 


^l2/K l 3 


1/31 


1/75 


1/33 



Tab. 7. 8. 2: Efficienza 



I valori riportati si riferiscono ad una analisi effettuata con tagli in ricostruzione rela- 
tivamente stretti sull'accettazione del vertice di decadimento della particella K s o K L . 
Cosi' facendo e' possibile ridurre il fondo dei decadimenti K L — > 7r°7r 7r° , a scapito co- 
munque deU'efficienza nel conteggio dei decadimenti K s — > 7r°7r e K L — > 7r°7r . 
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gamma 






3,4 


3,4 


2,5 


vertice 






200,1000 


200,5000 


-1000,5000 




in 


out 








K s -> 7r°7r° 








0.999 


0.923 


0.783 







1 


0.001 


0.077 


0.217 


in0/tot=0.666 


1 





0.820 


0.337 


0.098 


inl/tot=0.334 


1 


1 


0.180 


0.663 


0.902 




tot 


1 


0.061 


0.273 


0.446 


Kl -> 7r°7r° 








0.999 


0.908 


0.741 







1 


0.001 


0.092 


0.259 


in0/tot=0.603 


1 





0.921 


0.344 


0.104 


inl/tot=0.397 


1 


1 


0.079 


0.656 


0.896 




tot 


1 


0.032 


0.316 


0.512 


K L 7r°7r 7T 








0.000 


0.705 


0.659 







1 


0.005 


0.295 


0.341 


in0/tot=1.000 


1 





0.000 


0.000 


0.000 


inl/tot=0.000 


1 


1 


0.000 


0.000 


0.000 




tot 


1 


0.005 


0.295 


0.341 


K L 2/K L 3 


1/33 


1/196 


1/139 



Tab. 7. 8. 3: Efficienza 



I valori riportati si riferiscono ad eventi per cui in almeno una cella del calorimetro 
elettromagnetico venga rilasciata un'energia non nulla. 

Per i decadimenti Ks — > 7r°7r° occorre tenere conto che solo il 61% dei decadimenti generati 
arrivano sul calorimetro elettromagnetico senza essere bloccati dall'anticoincidenza Ks ■ 
Per i decadimenti Kl — > 7r°7r° e Kl — > 7r°7r°7r occorre considerare che gli eventi analizzati 
si riferiscono ad una piccola parte degli eventi effettivamente presenti nella realta', in 
quanto gia' a livello di generazione da parte del Montecarlo vengono creati soltanto eventi 
per cui il vertice di decadimento rientra in un intervallo di (—100, 3500 cm) dalla targhetta 
K s - 

Per il calcolo del doppio rapporto R da una successiva analisi off-line verranno presi in 
considerazione soltanto decadimenti K s — > n ir e K L — > 7r°7r° per i quali energia e vertice 
di decadimento rientrano entro determinati intervalli di accettazione. Cosi' facendo si 
puo' ovviare all'inconveniente di non conoscere esattamente la totalita' degli eventi di 
decadimento utili alia misura di R. 

Tali risultati suggeriscono, al fine di una buona efficienza nell'acquisizione dei decadi- 
menti K s — > 7r°7r° , K L — > 7r°7r° , senza una eccessiva accettanza del fondo K L — > 7r°7r 7r° , 
rispetto ai tagli proposti [8], un restringimento dell' intervallo di accettazione per il vertice 
di decadimento del K , e un allargamento dell'intervallo di accettazione per il numero di 
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gamma visti nella singola proiezione sul calorimetro elettromagnetico. 





proposta 


modifica 


gamma 
vertice 


(3,4) 
(200,5000) 


(2,5) 
(200,2000) 



Tab. 7.8.4: Tagli: proposta e modifica 



7.8.1 Variazione di soglia 

Una soglia di 1.0 GeV per il filtro temporale implica una diminuzione dell'emcienza totale 
nell'accettazione degli eventi buoni dello 0.53% per i modi di decadimento Kg — > 7r 7r e 
K L — > 7T°7r , mentre l'memcienza nell'accettazione dei decadimenti K L — > 7r 7r°7r° si 
riduce del 2.17%. 

7.8.2 Variazione di raggruppamento delle celle 

Variando il raggruppamento in blocchetti delle singole celle, e passando ad una configu- 
razione in cui ciascun blocchetto e' formato da 16x2 celle anziche' da 8 x 2 celle, nel caso 
in cui il valore della soglia del filtro temporale venga posto a 1.0 GeV, l'efncienza totale 
nell'accettazione degli eventi buoni cala del 1.37% per i modi di decadimento K s — > n°ir e 
K L — > 7r°7r° , mentre l'memcienza nell'accettazione dei decadimenti K L — > 7r 7r 7r° si 
riduce del 1.52%. 

7.9 Eventi accidentali 

Per una accurata simulazione della realta' fisica occorre prendere in considerazione, oltre 
alle fonti di rumore interne all'apparato, anche eventuali fenomeni che, dall'esterno, pos- 
sono disturbare o sporcare l'evento in esame. 

Tali fenomeni, comunemente indicati come eventi accidentali, costituiscono, insieme 
alia rimozione del fondo Kl — > 7r°7r°7r , il grosso problema a livello di una analisi dei dati 
a qualsiasi livello. 

Un decadimento ad esempio con soli tre gamma rivelati sul calorimetro elettromag- 
netico, in coincidenza del quale sia presente un gamma accidentale, puo' essere rivelato 
dall'apparato come un evento a quattro gamma, e, quantomeno a livello di numero di 
picchi ricostruiti, candidato ad essere un evento Kg — * 7r°7r° o Kl — > 7r°7r° , quindi da 
registrare. 
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Attraverso una analisi dei decadimenti K s — > 7r°7r° , K L — > 7r°7r° , K L — > 7r°7r 7r° in 
presenza di eventi accidentali, siano essi dovuti a gamma esterni particolarmente energetici 
non schermati dall'apparato, siano essi frutto di altri modi di decadimento del K, e' 
possibile verificare la variazione di efficienza del Trigger neutro rispetto ad una condizione 
ideale di funzionamento precedentemente trattata. 

Una analisi del comportamento del Trigger neutro in presenza di eventi esterni e' ef- 
fettuata con un rate di gamma accidentali pari al rate di decadimento di Ks e Kl che 
decadono in zona fiduciale, ed in modo tale da sincronizzare il gamma accidentale con la 
finestra del picco in energia rilasciata sul calorimetro formato dagli shapers. 

Per la distribuzione di tali eventi accidentali sono stati creati appositi moduli Fortran 
che, avvalendosi del metodo Montecarlo, generano un gamma con energia, impulso, punto 
d'impatto e tempo di impatto sul calorimetro elettromagnetico determinate 
Mentre per la distribuzione del punto di impatto sul calorimetro e' prevista una distri- 
buzione casuale uniforme su tutta la sua superficie, per quanto riguarda l'energia del 
gamma accidentale, questa, sperimentalmente misurata nel corso dei test effettuati su 
NA48, risulta distribuita con valor medio /x = 20 GeV, sigma a = 15 GeV, asimmetria 
con skewness pari a 2. 

L'interazione del gamma col Krypton liquido, e la conseguente produzione dello sciame 
elettromagnetico, vengono simulate, ai fini del calcolo dell'energia rilasciata su ciascuna 
cella del calorimetro, attraverso l'utilizzo di una distribuzione dell'energia in funzione della 
distanza dal centro dello sciame. Tale distribuzione, nota sperimentalmente e funzione 
dell'energia del gamma incidente, risulta dalla sovrapposizione di tre esponenziali decre- 
scenti in funzione della distanza e corrispondenti rispettivamente a testa, corpo, coda del 
singolo sciame elettromagnetico. 

II tempo di arrivo del gamma accidentale, in questo tipo di simulazione, viene scelto 
in modo tale da cadere, con distribuzione casuale uniforme, all'interno della finestra tem- 
porale del picco in energia rilasciata sul calorimetro elettromagnetico, e corrispondente a 
110 ns. 

Con la convenzione precedentemente adottata, ovvero indicando con inO e ml rispet- 
tivamente gli eventi in ingresso non buoni e buoni in assenza di accidentali, e con outO e 
outl gli eventi in uscita giudicati non buoni e buoni in presenza di accidentali, si otten- 
gono per i tre modi di decadimento i seguenti risultati: 



blocchetti 


8x2 


soglia 


0.5 


energia 


(60,180) 


baricentro 


(0,10) 



Tab. 7.9.1: Configurazione e tagli 
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gamma 






3,4 


2,5 


2,5 


vertice 






200,5000 


200,2000 


-1000,5000 




in 


out 








K s -> 7r°7r° 








0.916 


0.984 


0.886 







1 


0.084 


0.016 


0.114 


in0/tot=0.666 


1 





0.468 


0.743 


0.373 


inl/tot=0.334 


1 


1 


0.532 


0.257 


0.627 




tot 


1 


0.234 


0.096 


0.285 


Kl -> 7r°7r° 








0.893 


0.977 


0.807 







1 


0.107 


0.023 


0.193 


in0/tot=0.603 


1 





0.378 


0.830 


0.217 


inl/tot=0.397 


1 


1 


0.622 


0.170 


0.783 




tot 


1 


0.312 


0.081 


0.427 


K L 7r°7r 7T 








0.698 


0.977 


0.559 







1 


0.302 


0.023 


0.441 


in0/tot=1.000 


1 





0.000 


0.000 


0.000 


inl/tot=0.000 


1 


1 


0.000 


0.000 


0.000 




tot 


1 


0.302 


0.023 


0.441 



Tab. 7.9. 2: Efficienza in presenza di accidentali intensi 



In questo modo, dal raffronto dell'efficienza del Trigger neutro in tali condizioni di 
utilizzo con l'efficienza in assenza di accidentali, e dalla conoscenza del rate di decadimenti 
K s e K L e del rate di eventi accidentali, e' possibile stimare l'efficienza dell'apparato di 
trigger in condizioni reali di funzionamento. 

II rapporto outO/inl indica la percentuale di eventi buoni che, di un gamma ac- 

cidental, non vengono accettati dal trigger. II rapporto outl/inO evidenzia la percentuale 
di eventi non buoni in ingresso che vengono invece giudicati buoni dal trigger in presenza 
di un evento accidentale; uno dei possibili casi in cui si puo' verificare tale situazione e', 
ad esempio, un decadimento in 7r°7r° con soli tre gamma sul calorimetro elettromagnetico 
(mO) piu' un quarto gamma accidentale che apporta una correzione alia ricostruzione di 
energia, baricentro dell'energia e vertice di decadimento, altrimenti scartati, tale da far 
accettare l'evento come decadimento in 7r°7r° . 

Detto r d il rate di decadimenti K s e K L , e r a il rate di eventi accidentali sul calorimetro 
elettromagnetico, la percentuale di eventi buoni contaminati da un evento accidentale e' 
data da: 

[7.10] P = -{2rr d -r a )=2rr a 

dove r e' la larghezza temporale del segnale di shaper. 
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Dalla conoscenza quindi dell'efficienza eo in assenza di eventi esterni, dell'efficienza in 
presenza di eventi accidentali come precedentemente definita, l'efficienza in condizioni 
reali di funzionamento e del Trigger neutro nell'acquisizione di eventi buoni risulta essere: 

[7.11] e = e -pio+Poi 

dove pio e' la percentuale di eventi buoni sporcati da eventi accidentali e non accettati dal 
trigger, e p i e' la percentuale di eventi non buoni resi accettabili dal gamma accidentale. 
I valori di p w e p i risultano dati da: 

7.12 pi0 = P . 

ml 

dove i rapporti relativi alle inefficienza dell'apparato si riferiscono ai valori calcolati in 
presenza di accidentali con rate equivalente al rate degli eventi di decadimento. 

In condizioni di rate di decadimenti di K$ e Kl pari a 1 MHz, e rate di eventi acci- 
dentali pari a 10 KHz, la percentuale di eventi contaminati da un gamma accidentale e' 
pari a 

[7.14] P = 2.2 • 1(T 3 

ovvero la probabilita' di contaminazione di un evento da parte di un accidentale e' tra- 
scurabile ai fini dell'esperimento. 

Applicando i tagli sulle quantita' fisiche ricostruite precedentemente introdotti (tagli 
modificati), si ottiene una efficienza di funzionamento dell'apparato di trigger, calcolata 
come rapporto tra gli eventi giudicati buoni in uscita e gli eventi effettivamente buoni in 
ingresso: 

K s -> 7r°7r° : e = 0.754 

K L -> 7r°7r° : e = 0.384 
Questo parametro indica l'effettiva efficienza dell'apparato di trigger in presenza di tali 
tagli, senza tenere conto della possibile dispersione dei gamma, e quindi del fatto che 
molti eventi giungono con un numero inferiore di gamma rispetto a quelli prodotti sul 
calorimetro. 

Per quanto riguarda invece i decadimenti di tipo K L — > 7r°7r 7r° , una percentuale di 
eventi pari a 6.09% viene candidata dal trigger neutro ad essere frutto di un decadimento 
in 7r°7r° . 
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